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tionship for tensile specimens extended at temperatures ranging from 196° to 
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Specimens of ingot iron in various initial conditions were extended to a specified 


strain at a selected temperature and subsequently extended to fracture at a different tempera- 


ture. 
corresponding values for a single-stage 
prestrain of the specimen increases 
The rheotropic embrittlement”’ 
aging 
at or below TRC 
in tension at room temperature 
room temperature 


test 


and 


decreases the ductility 


work 
vary with the heat treatment and prestrain history of the ingot iron 


The deviation of the true stress-strain values for the second stage of these tests from 
generally 


increase at a decreasing rate as the 


hardening (strain hardening plus strain 


The brittleness 


of normalized or hot-rolled ingot iron is partially curable by prestraining 
However, prestraining in tension of annealed ingot iron at 
retained at 154 


and 196°C The predomi 


nance of strain aging at the higher temperatures and strain hardening at the lower tempera- 
tures causes the total work hardening of the specimen during deformation to maximum load 


to increase as the 


l. Introduction 


In recent vears several investigations |1 to 11]? 
on the mechanical properties of metals at low tem- 
peratures have shown conclusively that the true 
stress-strain relationship of metals and alloys ex- 
tended in tension depends on the strain-temperature 
history, as well as on the instantaneous values of 
strain, strain rate, and test temperature. In most of 
these investigations very little information was re- 
ported on the part played by strain aging during the 
prior history of the specimen, as strain aging was not 
pronounced in most of the metals studied. A previ- 
ous investigation [10] at the Bureau showed that 
strain aging affected the true stress-strain relation- 
ship of ingot iron specimens extended in tension at 
slightly elevated temperatures and even at some sub- 
zero temperatures 120° to 100° C). Within the 
temperature range, —196° to +-100° C, of that in- 
vestigation the strain hardening of ingot iron in- 
creased with a decrease in temperature, whereas the 
strain aging decreased with a decrease in temperature. 
Therefore, work hardening * of ingot iron may either 
increase or decrease with increase in testing tempera- 
ture, depending upon the relative magnitude of the 
changes in strain hardening and strain aging with 
Variation in temperature. 

The rate of strain hardening of ingot iron is in- 
fluenced by the mechanism of deformation. The 
pacing of the slip zones within the crystals and also 
the number of slip bands in a zone vary with the tem- 
perature of deformation of the metal. These factors 
are believed to be directly associated with the increase 
n the rate of strain hardening of the metal with de- 
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Army, through the Ordnance Corps, Watertown Arsenal 
cures in brackets indicate the literature references at the end of this paper 
As used in this paper the term “‘work hardening” relates to the increase in 
tress with true strain and is the combined effect of strain hardening and 
wing lhe term “strain hardening” refers to that portion of the hardening 
rengthening of the metal during plastic deformation that is not time depend- 
' and occurs concurrently with strain, whereas the term “‘strain aging’’ refers 
hat portion of the hardening or strengthening occurring during or subsequent 
istic deformation that is time dependent and occurs as a result of previous 
ind not concurrent strain 
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test temperature is either increased or decreased from about 120° ¢ 


crease in testing temperature Moreover, as pre- 
viously reported [10], considerable deformation by 
twinning was observed in specimens of ingot iron 
slowly extended in tension at temperatures below 
about —120° C. The extent of the twinning during 
the tension test increased greatly as the temperature 
was lowered from about 120° to —196° C. The 
influence of this factor on the rate of strain hardening 
of the iron may be considerable 

The purpose of the present investigation was to 
obtain additional information on the influence of the 
strain-temperature history, and especially the effect 
of strain aging, on some of the tensile properties and 
the true stress-strain relationship for specimens of 
ingot iron subsequently tested in tension. 


2. Material 


The material used in this investigation was ingot 
iron in the conditions as annealed, hot-rolled, 
quenched and tempered, normalized, and cold-drawn. 
A detailed description of this material was reported in 
a previous paper [10]. The principal chemical con- 
stituents (in percent by weight) other than iron are 
as follows: Carbon, 0.02; manganese, 0.02; phos- 
phorus, 0.005; sulfur, 0.018; silicon, 0.002; copper, 
0.10; oxygen, 0.058; nitrogen, 0.002; and hydrogen, 
0.0005. This iron was prepared from a single melt and 
was furnished by the manufacturer in the conditions 
as hot-rolled and as cold-drawn to 14- and 24-percent 
reductions of area. The heat treatments applied to 
the ingot iron in this laboratory are summarized in 
table 1. 

Cylindrical tensile specimens with a 2-in. gage 
length were used. The reduced section was gradually 
tapered from each end; the diameter (0.438 in.) at 
the midsection of the gage length was about 0.003 in 
than at at each end. The specimens were 
finished to the final dimensions by grinding and 
polishing in the axial direction. The ends of the 


less 


| specimen were machined with % in. X10 threads, and 
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the shoulder fillets were machined to a radius of 0.75 
in 
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3. Method of Investigation 


The specimens were extended in tension in a pen- 
dulum hydraulic testing machine of 50,000-lb 
capacity. The specimens, except tested at 
room temperature, were fully immersed in an appro- 
priate liquid maintained at the desired temperature 
The loading was so controlled that the rate of reduc- 
tion of area of the specimen, bevond the region of 
initial yielding, was maintained at approximately 1 
percent/minute. Measurements of the diameter ot 
the specimen throughout the course of each test were 
made by means of a specially designed reduction of 
area gage [12]. Changes in diameter of a specimen 
were measured by this instrument with an accuracy 
of +0.0001 in. A detailed description of the testing 
equipment and the method of maintaining the desired 
temperature is given in previous papers [10, 12]. 


those 


PaBLE 2 Two-stage test conditions 


remperature of specimet 


Load 

bet ween 

First stag Second stage Betw stages 

rv 
( ‘ ‘ , 

1 7s s a” 
Roon 7s s 200 
Deo 1m | 200 
Deo iM ha | 200 
I> lus lw a 
7s Room - 0 
!* Ive ” a0 
1w Room 1m 0 
1M do 15 0 
1” do 10 0 
ive 7s 7s 200 


i Two special tests also were made in which the load obtained at the end of the 
first stage was maintained between stages 
? Specimen was removed from the testing machine bet ween stages 


In the present investigation, a specimen was 
extended at a selected temperature to a specified 
strain value and subsequently extended to fracture at 
a different temperature. In general, the specimen 
was kept in the testing machine between stages under 
mM email load of approximately 200 lb to maintain the 
alinement of the specimen in the testing machine. 
However, when the second stage of a two-stage test 
was at room temperature, the specimen was removed 
from the testing machine between stages so that the 


adapters and adjacent portions of the head and ba 
of the testing machine could be quickly brought 
room temperature. Moreover, in two special t 
the load attained at the end of the first stage w 
maintained on the specimen in the machine betwe: 
stages The specimens were kept cool betwe 
stages, as reported in table 2, in order to minim 
strain aging during this period. Usually the seco 
stage was begun approximately } hr after complet 
of the first stage. 


4. Results and Discussion 


True stress-strain data obtained in two-stage t 
sion tests on ingot-iron specimens at selected temp: ' 
atures ranging from 196° to 100° C are su 
marized in figures 1 and 2, and 4 to 9. True stress- 
strain curves obtained in the single-stage tension 
tests, reported in a previous paper [10], are includ 
in these graphs as broken-line curves, and the upper 
end of these curves represents initial fracture. Th 
true stress-strain values representing the completion 
of the first stage and the data obtained in the second 
stage of the two-stage tension tests are shown. Thy 
data obtained in the first stage of the test conforms 
to the curve representing the corresponding deforma- 
tion in a single-stage test and these values are, in 
most cases, not plotted in these graphs. The tru 
stress-strain values during the unloading and th: 
initial part of the reloading also are not shown in som: 
of the figures. The true-stress and true-strain values 
attained at the fracture area of the specimen at th: 
initiation of fracture are designated in this paper as 
the fracture stresses and fracture strains, respectivel) 
and, except for certain specifically designated cases 
they are represented by the final true stress-strain 
values plotted in figures 1 and 2, and 4 to 9 


4.1. Annealed Ingot Iron 


In order to determine the direct effect of a reduc- 
tion of load in a tension test for a short time on the 
subsequent true stress-strain values, a special tension 
test was made at a single temperature (room temper- 
ature) on a specimen of annealed ingot iron. The 
deformation of this specimen was carried out in four 
stages and included an initial loading to a selected 
extension, an unloading to 200 lb for a period of 5 
minutes, followed by a reloading to a second selected 
extension, then repeating the procedure until a total 
of four stages had been completed. As shown in 
figure 1, the true stress-strain values for this four- 
stage test generally coincided with those of a single 
stage test made on a specimen of annealed ingot iron 
at room temperature, except for the true stresses at 
the initial strains of the second, third, and fourth 
stages; these true stresses were slightly higher than 
the corresponding stresses in a single-stage tes! 
Apparently, the combined effect of the strain aging at 
room temperature during the 5-minute periods 
between stages and the unloading and reloadii 
was sufficient to increase the true-stress values at 1! 
initial yielding in the last three stages but was n 
sufficient to affect appreciably the subsequent tr 
stress values. It is believed that, in general, sh: 
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strain values obtained with an annealed 


es 1. True strese- 
/ specimen that was extended in a four-stage tension 


fem perature 


interruptions of this type in a tension test with ingot 
iron at room or lower temperatures, in which the 
temperature of the specimen is held constant during 
the test. have very little effect on the true stress- 
strain values, except those during the initial yielding 
n the next stage of the test 
Four specimens of annealed ingot iron were ex- 
ended in tension at room temperature to specified 
strains (indicated as short vertical bars on the 27° C 
vein fig. 2), and subsequently extended to fracture 
it 196° C. As shown in figure 2, these specimens 
exhibited very little ductility in the second-stage 
test, although a specimen of annealed ingot iron that 
was extended in a single-stage test at 196° C 
deformed to a strain of about 0.05 before initial frac- 
In the two tests, in which the deformation 
of the specimen in the first stage of the test at room 
temperature was greater than the maximum load 
strain, the specimens fractured in the second-stage 
test without any appreciable plastic deformation, 
and the fractures occurred at the shoulder of the 
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at room temperature and — 196° C 








Portions of two specimens of annealed ingot tron 


FicuRE 3 
that were extended in tension at room te mperature to strains 


beyond marimum load and subse quently Jract ired in tension 
at 196” ¢ 


A, True strain of 0.62 at room temperature; B, true strain of 0.94 


necked section (fig. 3 The final true stress-strain 
values as determined by the final load and minimum 
cross-sectional-area measurements for these speci- 
mens are designated as points A and B in figure 2. 
As the fractures did not occur at the minimum 
cross section of the specimens, the stress values 
represented by points A and B are not fracture 
stresses; these data indicate only that the fracture 
stresses at — 196° C for the position of the minimum 
section of the specimens of annealed ingot iron that 
were extended in tension at room temperature to 
these designated strain values were greater than the 
indicated stress values The fracture stresses and 
fracture strains at —196° C of these specimens under 
the multiaxial stress conditions existing at the posi- 
tion of the actual fracture in the shoulder of the 
necked section are shown by the corresponding 
points A’ and B’ in figure 2. It may be noted that 
the fracture stresses at 196° C in the shoulder 
portion of the neck of the specimens are lowered to 
values below that of the tension specimen fractured 
at 196° C in uniaxial tension after extension at 
room temperature to maximum load. This may be 
attributed to the influence of biaxial compressive 
stresses at these positions. 

As the ductility of specimens of annealed ingot 
iron in single stage tests at 196° C was small 
(strain of approximately 0.05), very little information 
could be obtained on the effect of prior strain at 

196° C of specimens of annealed ingot iron on 
their tensile properties at room temperature. The 
data obtained on a specimen extended in the first 
stage of the test at —196° C to a strain of about 0.02 
and then extended to fracture in the second-stage 
test at room temperature indicated that the total 
work hardening during the deformation at 196° C 
was approximately equivalent to the total work 
hardening of a specimen deformed the same amount 
at room temperature; the true stress-strain values 
obtained in the second stage of the test approxi- 
mately coincided with the corresponding values ob- 
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tained in the single-stage test at room temperature 
2 Thus the increase in the strain hardening 
above that at temperature was 

equal to the strain room 


(fig 
ut 196° © room 
approximately aging at 
temperature 

\ specimen of annealed ingot iron was prestrained 
temperature to maximum load 
to fracture at 154° C. This 
specimen retained considerably ductility and 
hence had a lower fracture stress specimen 
extended to fracture in a single-stage test at 154°C 
The strain during the second 


mn tension at room 


ind then extended 
less 


than a 


as shown in figure 4 
stage of this test was very small in comparison with 
the strain beyond maximum load in the single-stage 
test at —154° C. The strain aging of the specimen 
during the deformation at temperature ap- 
parently made the specimen brittle at the 
lower temperature 


room 
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with unneated 
fwo-stage 


ss-strain values obtained 


extended in 


154° ¢ 


Fiaure 4 True stre 


specimens that were tension 


ingot-iron 

tests at room temperature and 

An annealed ingot-iron specimen was prestrained 
in tension at 154° C to maximum load and then 
extended to fracture at room temperature. The 
true stress-strain curve for the second stage of this 
test (fig. 4) is considerably below the corresponding 
part of the curve for a single-stage test with a speci 
men at room temperature, and thus the total work 
hardening of the specimen during the first stage of 
the test at —154° C was less than that of a specimen 
deformed to the same strain at room temperature 
Strain aging in specimens of ingot iron during defor- 
mation in tension at 154° © is very slight [10] 
Moreover, strain hardening increases with a lower- 
ing of the test temperature from room temperature 
to —154°C. ‘The data in figure 4, therefore, indicate 
that the strain aging of the iron at room temperature 
was much greater than was the increase in strain 


from 


hardening with decrease room tempera 
to 154° C 

In order to determine the specific influence of 
degree of the prestraining of annealed ingot- 
specimens at one temperature on their deformat ) 
characteristics at another temperature, tWo seri 
One series of specim 


two-stage tests were made 
was prestrained In tension at room temperature 


different strains and subsequently extended 
fracture at 78°C. The true stress-strain data 
these tests are summarized in figure 5. The t 


stress-strain curves for the second stage of thes 
tests are above the corresponding portion of the cu 
Another series of 


7a” ¢ 


for a single-stage test at —78° C 
specimens Was prestrained in tension at 
to different strains and subsequently extended 








fracture at room temperature The true stress 
strain curves (fig. 5) for the second stage of each 
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Figure 5 True stress-strain values obtained with ann 
ingot-tron specimens that were extended in two-stage ter 
tests at room temperature and r8° ( 


of these tests, except for the specimen prestrained 


only to the lower yield, are below the corresponding 


portion of the curve for a single-stage test at room 
temperature. The curve for the second-stage tes 
of the specimen prestrained to the lower vield 
coincided with the curve for the single-stage test at 
room temperature. The combined effect of the strain 
hardening and strain aging during the prestraining at 
78° C in this test was equivalent to the combined 
effect of the strain hardening and strain aging 
during the deformation to the same strain at roo! 
temperature. In general, the deviation of the tru 
stress-strain curve for the second stage of a two 
stage test from the curve for a single-stage tes! 
the same temperature (fig. 5) increases at a «i 
creasing rate as the amount of prestraining in |! 
The above data indicate (hat 


first stage increases 
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the strain aging of a specimen of annealed ingot 
ron at any strain during deformation in tension at 
a constant temperature with the strain; 
the straim aging at any strain depends greatly upon 
the amount of the prior strain aging of the specimen 

In order to determine the effect of maintaining a 
relatively high stress on the specimen between test 
stages, two special tests were made in which the load 
attained at the end of the first stage was maintained 
specimen until the beginning of the second 
stage The true values obtained in 
these tests (fig. 5 
the true stress-strain curves of corresponding two- 
tests in which the load was reduced to about 
200 Ib between stages. The retained ductility at 

78 C of the specimen prestrained at room 
temperature to a strain of about 0.26 with no re- 
duction of the load between slightly 
greater than that of the specimen prestrained at 
temperature to approximately the same 
with the load reduced to 200 lb between 
This small difference in ductility, however, 


Varies 


on the 
stress-strain 
show no significant deviation from 


stages Was 


room 


amount 


stages 
is not believed to be significant In general, the 
strains to fracture of the specimens in the two- 


stage tests summarized in figure 5 were approximate- 
ly the same at both temperatures and also were 
about the same as those for the single-stage tests at 
these temperatures 

Ripling and Baldwin in reporting on the phe- 
nomenon they called the ‘“rheotropic embrittlement” 
of steel [11] concluded that a large part of the de- 
ficiency in ductility of annealed steels at low temper- 
atures is rheotropic and can be removed by prestrain- 
ing under more ductile conditions at a higher temper- 
ature. The data obtained in the present inves- 
tigation on annealed ingot iron, however, indicate 
that the low-temperature embrittlement of this 
material is not rheotropic in the manner proposed 
by Ripling and Baldwin; prestraining of the annealed 
ingot iron at room temperature reduced, rather 
than increased, the ductility retained either at 

154° or 196° C. The embrittlement due to 
strain aging during the prestraining of the specimens 
at room temperature, is believed to be the main 
cause for the loss in the retained ductility at —154 
ol 196° ©. Another factor that probably has 
considerable influence on the retained ductility is 
a decrease in the degree of deformation twinning 
occurring at the low temperature following a pre- 
straining at room temperature. This latter factor 
is being studied in some detail and will be discussed 
in #& separate report 


4.2. Normalized Ingot Iron 


The true stress-strain data for two normalized 
ingot-iron specimens that were extended in tension 
" two-stage tests at room temperature and —78° C 
are presented in figure 6. The effect of the prestrain- 
ing in the first stage of the tension test on the total 
work hardening of these specimens was similar to 
that already described for two-stage tests with an- 
nealed ingot iron at these temperatures. The duc- 
tility, however, was affected quite differently. \s 








Frat RE 6 True stress-strain va s obtained with normalized 


ingot-iron specimens that were extended in two stage tension 
tests at room temperature and gs” ¢ 


illustrated in figures 25 and 26 of the previous paper 
[10], the ductility of normalized ingot iron, and to a 
lesser degree that of hot-rolled ingot iron, extended to 
fracture in single-stage tests was much lower at 

78° C than at room temperature. The strain to 
fracture of ingot iron in the annealed, quenched, and 
tempered, or cold-drawn conditions, however, was 
nearly the same at room temperature and at —78° C 
The decrease in the ductility at —78° C of the ingot 
iron in the normalized or hot-rolled condition may 
be associated with a greater amount of carbon and 
nitrogen initially retained in solid solution in these 
conditions than in the other conditions. The data 
for the two-stage tests of the normalized ingot iron 
(fig. 6) show that the iron in this condition apparently 
is subject to rheotropic embrittlement in the manner 
proposed by Ripling and Baldwin; prestraining the 
specimen at room temperature to just beyond the 
maximum load increased greatly the ductility re- 
tained at —78° C, as indicated by the fracture strain 
of approximately 1.1 for this two-stage test and the 
fracture strain of approximately 0.5 for the single 
stage test at —78° C.  Prestraining a specimen at 

78°C to the maximum load decreased the ductility 
retained at room temperature. Apparently, the 
rheotropic embrittlement of the specimen during the 
prestraining at —78° C was not completely removed 
by the subsequent deformation at room temperature 


4.3. Quenched and Tempered Ingot Iron 


The data obtained in tWo-stage tension tests with 
specimens of quenched and tempered ingot iron are 
summarized in figure 7. These results indicate that 
the effect of preStraining specimens at room tempera- 
ture and —78° C is approximately the same as that 
described previously for specimens of the annealed 
ingot iron; the difference in magnitude of strain aging 
at —78° C and room temperature was greater than 
the difference in strain hardening at these tempera- 
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Fiovurnre 7 True etreas-strain va s obtained th quen hed and 
tempered ingot on specimens that were extended in two-stage 

fenasion teata at oom temperature ind y f 


tures, and the predominating effect of the strain 
aging resulted in greater total work hardening during 
the prestraining at room temperature than at 78°C 
The ductility retained by the specimens in the second 
stage of the tests was not affected appre iably by the 
difference in the prior strain histories 


4.4. Hot-Rolled Ingot Iron 


As the ductility of the hot-rolled ingot iron was 
considerable, even at low temperatures (fracture 
strain of approximately 0.3 at 196° ©), specimens 
were prestrained in tension at temperatures ranging 
from 196° to +100° C-and then extended to frac- 
ture at different temperatures. The true 
strain data for these tests and the curves obtained in 
single-stage tests at these temperatures are summar- 
ized in figure 8. Prestraining specimens approxi- 
mately to maximum load at 196 120", or 

78° C had little specific effect on the ductility 
retained in the second stage of the tests at room 
temperature. The deformation under ductile condi- 
tions at room temperature apparently removed any 
rheotropic embrittlement that may have occurred 
during the prestraining at low temperature. As the 
work hardening of a specimen of ingot iron during 
deformation at —196° C is mainly strain hardening, 
and as strain hardening increases with decrease in 
temperature from —78° to — 196° C, the data showing 
equal work hardening of the hot-rolled ingot-iron 
specimens during the prestraining at 78° and 

196° C (fig. 8) indicate that strain aging occurred 
to a considerable extent during the prestrain of the 
specimen at —78° C. The work hardening during 
the prestraining in tension at — 120° C of a specimen 
approximately to maximum load was the least for 
this series of tests and resulted in the greatest lower- 
ing of the points representing the true stress-strain 
values for the second stage of the tests from the curve 
for the single-stage test at room temperature. The 


stress- 


strain aging of hot-rolled ingot iron during defo 
tion was very much less at —120° C than at —7s ( 
However, as shown in figure 12 of a previous in, 
gation [10], strain aging occurs to an appre: 
extent in specimens of hot-rolled ingot iron defor: od 
in tension at 120° C and at 138° C. 

Although the total work hardening (mainly strvi) 
hardening) of a specimen during prestraining «p 
proximately to maximum load at — 196° C was ab 
equal to the total work hardening (strain hardening 
plus strain aging) of a specimen deformed to the sam, 
strain at —78° C, the resulting rheotropic properties 
of these two specimens were different. The points 


representing the true stress-strain values for thy 
second-stage test at 78° C of the specimen pre 


strained at —196° C gradually rise above the corr 
sponding portion of the curve for the specimen ex 
tended in the single-stage test at 78° C (fig. 8 
These data indicate that the strain aging during th. 
deformation from maximum load to fracture of th 
former specimen was much greater than that of the 
latter specimen. As the fracture strains of the abov: 
specimens were about equal, the deformation during 
the second-stage test at —78° C apparently removed 
most of the low-temperature embrittlement that had 
occurred during the prestraining at 196° C 

The specimen of hot-rolled ingot iron that was pre- 
strained in tension at room temperature to a small 
strain of approximately 0.03, broke in the threads 
during the second-stage test at 196° C without any 
appreciable plastic deformation of the specimen 
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The final value plotted in figure 8 for this test does 
not represent the fracture stress and fracture strain 
at 196° C of the prestrained metal. As the yield 
stress of this specimen at — 196° C was at least equal 
to or greater than that represented by the final 
plotted value, and as this point lies approximately on 
the true stress-strain curve for the single-stage ten- 
sion test on hot-rolled ingot iron at 196° C, the 
prestraining of the specimen at room temperature, 
apparently did not decrease the stress necessary for 
continued deformation of the metal at 196° C 

Specimens of hot-rolled ingot iron that were pre- 
strained in tension at room temperature approxi- 
mately to maximum load and then extended to 
fracture at 78”, 120°, and 196° C, exhibited 
strain-aging effects similar to those described above 
for the specimens prestrained at low temperatures. 
However, a difference was observed in the effect of 
the prestraining at room temperature on the ductility 
retained at lower temperatures, and this should be 
pointed out. Prestraining at room temperature 
approximately to the maximum load improved the 
ductility retained by the specimen at the lower tem- 
perature and the fracture strains of the specimens 
in the two-stage tests were greater than those of the 
specimens fractured in the single-stage tests at the 
same low temperatures fig. 8) Thus, the rheo- 
tropic embrittlement of the hot-rolled ingot iron 
during the second stage of these tests was reduced by 
the prestraining under the more ductile conditions at 
room temperature 

\ specimen of hot-rolled ingot iron was prestrained 
in tension at 100° C to a strain slightly beyond the 
maximum load and then extended to fracture at 

78°C. The strain aging of this specimen during 
the deformation at 100° C was much greater than 
that of specimens deformed at room temperature, 
and is indicated by the position of the true stress- 
strain values in figure 8 for the second stage of the 
test; the elevation of these points above the curve for 
the single-stage test at —78° C is much greater than 
that for the specimen prestrained at room tempera- 
ture. The fracture strain of the specimen prestrain- 
ed at 100° C was about equal to that of the specimen 
extended to fracture in a single-stage test at —78° C. 
However, it was considerably less than that of the 
specimen prestrained at room temperature and ex- 
tended to fracture at 78° CC. The embrittlement 
of the specimen by the great amount of strain aging 
during the deformation at 100° C apparently was 
about equal to the rheotropic embrittlement during 
the deformation of a specimen to the same strain in 
the single-stage test at —78° C, 

The true stress-strain data for hot-rolled inget 
iron specimens (fig. 8) indicate that the strain hard- 
ening is the predominating factor affecting the var- 
‘ation of the total work hardening with temperature 
at temperatures below —120° C, whereas the strain 
aging is the predominating factor affecting the total 
work hardening at temperatures above —120° C. 

As described earlier, the ductility of specimens of 
normalized ingot iron was much less at —78° C and 
at lower temperatures than at room temperature, and 


prestraining of a specimen at room temperature 
greatly improved the ductility retained at —78° C 
The ductility of the annealed ingot. iron specimens 
was less at 154° C and —196° C than at room 
temperature. However, prestraining specimens of 
annealed ingot iron at room temperature reduced 
the ductility retained at 154° C and 196° C 
Moreover, the ductility of specimens of ingot iron 
in the condition either as annealed, or as quenched 
and tempered, was about the same at —78° C as at 
room temperature and prestraining specimens at 
room temperature had little, if any, effect on the 
ductility retained at 78° C The effect of pre- 
straining of ingot iron at room temperature on the 
rheotropic embrittlement of the specimens at lower 
temperatures apparently depends greatly upon the 
initial condition of the iron 


4.5. Cold-Drawn Ingot Iron 


The true stress-strain data obtained in two-stage 
tension tests with specimens of cold-drawn ingot iron 
are summarized in figure 9. These samples of the 
originally hot-rolled ingot iron had been cold-drawn 
by the manufacturer to 14- and 24-percent reduction 
of area, and were aged at room temperature for 
several months. The total work hardening during 
the prestraining in tension at room temperature to 
a strain of approximately 0.25 of a specimen of ingot 
iron cold-drawn to 14 percent reduction of area was 
slightly greater than that of a specimen deformed to 
the same strain at —78°C. The total work harden- 
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True stress-strain. values obtained with*cold-drawn 
ingot-iron specimens that were extended in two-stage tension 
tests at room temperature and 78° ( 
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ing during the prestraining in tension to a strain 
~ 


of approximately 0.25 at room temperature and at 


78° C of specimens of ingot iron cold-drawn to 
24-percent reduction of area was approximately the 
same as indicated by the coimecidence of the true 
stress-strain values for the second stage olf these 
tests with the true stress-strain curves for the single 


of strain aging during the 
specimens of the 


The degree 
tension of 


tage tests 
cdetormation in these 
cold-drawn iron was much less than that of the hot- 
olled ingot-iron specimens due to the prior strain 


that had occurred during t he cold-drawing and 


nein 
storage of the tron 
5. Summary 
\ study was made to determine the effect of the 


strain-temperature history of ingot iron on the true 
stress-strain relationship for specimens extended in 
from 196” to 

Loo U Specimens of annealed normalized, 
quenched and tempered hot-relled and cold-drawn 
ingot iron were extended in tension to specified true- 
strain values at 


quently extended to fracture at 


tension at temperatures ranging 


a selected temperature and subse- 
a different temper- 


ature. The true stress-strain values obtained in 
these two-stage tests were compared with those 
obtained in single-stage tension tests at the same 


the second stage of these tension 


temperatures deviation of the true stress- 
strain curves for 
tests with ingot iron from the true stress-strain curve 
for a single-stage test at the same temperature gen- 
erally increases at a decreasing rate as the prestrain 
of the specimen is increased 

The amount of strain aging and the rheotropic 
brittleness of specimens of ingot iron vary with the 
heat treatment and prestrain history of the speci- 
mens. The rheotropic brittleness at —78° C or at 
lower temperatures of normalized or hot-rolled ingot 
iron Is sensitive to and in part curable by prestrain 
under ductile conditions at higher temperatures; the 
ductility retained by the specimens in subsequent 
extension in tension at low temperature is increased 
However, prestraining specimens of annealed ingot 
iron in tension under ductile conditions at room 
temperature decreases the ductility retained in sub- 
sequent extension in tension at 154” o1 196° C 

The work hardening of ingot iron during deforma- 
tion in tension at temperatures below 120° © is 
mainly strain hardening, and increases with lowering 
of the testing temperature. As the temperature is 
increased above — 120° C, the influence of the strain 
aging during the deformation in increasing the total 
work-hardening rate predominating 
factor, and the rate of work hardening increases with 
increase in testing temperature. The data obtained 
on the tension specimens of hot-rolled ingot iron 
indicate that the part of the total work hardening 
due solely to the strain aging during the deformation 


becomes the 
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to maximum 
as low as —78° C 


load is considerable at temperat 
Even at —78° C the work har 
ing due to the strain aging is approximately equ 
the increase in strain hardening that occurs as 
testing temperature is lowered from —78° to — 196 
At 100° C, the magnitude of the strain aging 
at least equal to or greater than the strain harden 
The data presented in this paper, together w 
those reported previously !10], show that the st: 
aging of ingot iron is a very important factor affect 
the true stress-strain relationship at temperatu 
as low as 120°C. Strain aging apparently oc 
to some extent at temperatures as low as —154° ( 
These data are contrary to the generally prevail 
belief that the strain aging of 
nificant at sub-zero temperatures 


ingot iron Is ins 
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Gradient Methods in the Solution of 
Systems of Linear Equations 


Marvin L. Stein 


rhe results of various experiments with iterative methods for solving systems of lineat 
ilgebraic equatio are discussed Modifications of the optimum grad ( <i are ¢ " 
vmared, and the rather interesting self acceleration properties of a class of ethods here named 
almost optir radient methods are pointed out 


l. Introduction 


The method of steepest descent, or the optimum 
radient method, has been known to mathematicians 
since the tume of Cauchy |1 Others who have dis- 
sed this method include Curry |2], Forsythe and 
\lo n [3], Householder |4], Kantoroviteh [5], and 
lemple (6). Its infrequent application in compu- 
tiona! work is no doubt due to the slowness with 
vhich it converges This slowness of convergence 
s unfortunately generally true of gradient methods 
However, with the advent of large-scale computing 
lachinery it has become feasible to seriously considei 
em in practi al numerical analvsis 
In a forthcoming paper, Hestenes and Stein {7 
SCLISS laree class of eradient proc edures for soly ine 
ems ol lineal equations The se procedures con- 
ain the optimum gradient method as a special case 
The present note is mainly a report on some nume: 
il experiments with them that were carried out on 
IBM Card-Provrammed Electronic Calculator at 
the Institute for Numerical Analysis of the National 
Bureau of Standards 
an experiment in which the problem of solving a 


Some attention is also rivel 
system o! linear equations was ¢ hanged to an equi 
pent ¢ renvalue problem and then solved bv a mod 
fication of one of the gradient methods discussed by 
Hestenes and Karush [8 The most striking result 
ne ¢ xperuments indicates that there is a large class 
eradient methods that is se/f-a ‘celerating, that is, 
which irregularly shows a large increase in the rate of 
mvergence without the introduction of any modi 
ication in the computational routine. This be- 
havior is in sharp contrast to that of the method of 
steepest descent, which in the light of the present 
results can no longer be considered as optimum from 
in over-all point of view unless modified by some 
special accelerating routine |9 


2. Summary of the Theory 


Let A denote an m-rowed and n-columned matrix, 
ind let 6 and « be m-rowed and n-rowed column 
ctors, respectively The 
then an m-rowed column vector 


vector rir) b Ay Is 
In the following 
tar , affixed to the symbol for any vector or 
atrix will indicate its conjugate transpose. Hence, 


if JT isa positive m m Hermitian matrix, we have 
H*=H r*Hr>0 


provided r#0. Clearly, the 
nonnegative real-valued function 


f(a "Th 1) 
vanishes if and onlv if 
iz—h, 2) 


Therefore, the problem of minimizing (1) is equiva- 
lent to the problem of solving the svstem (2), pro 
vided a 


vector 7 


solution exists. If no solution exists, a 
which minimizes f(r) vields a best fit of 
b by Ax in a least-squares sense with the metrir 
determined by H 

Hestenes and Stein 7| have analyzed the following 
algorithm for constructing the minimum of (1 
Consider iterations ot the ty pe 


J I4 ar, e ) i . . , > 


is chosen initially and where, after z,; has 
been determined the gra lient vector & is defined by 
the rule A*®Hr(a If ¢ 0, the problem is 
solved If =. +0, the sealar a, is taken to be of the 
form a 3,9 


W here J 


whe re 
A*HTA$ 
and 8; is any complex number. The 


has been shown to converge to the minimum of 
fia provided the coefficients 8, satisfv the conditions 


sequence 3 


where 6 is arbitrary on the range 0<6< 1. 

Setting 8 | gives the optimum gradient method 
which has the following geometrical interpretation 
Starting at z, one proceeds along the normal to that 
member of the family ot concentric ellipsoids 


Ae constant, } 


which passes through z, until a point tangent to 
another ellipseid of the family is reached. One thet 
goes along the norma! to this secend ellipsoid unti! 
one is again at a point of tangency to a member of (4 
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and so on until the common center is reached. Clear- 
lv, it would be highly desirable to land on an axis 
of the family. However, as simple two-dimensional 
examples illustrate, one usually overshoots the major 
axis by proceeding all the way to a point of tangency 
Hence Hestenes that 
of the allowed values of 8, some of chosen 
from the range 8 l would yield better convergence 
than 8,=1. This conjecture is the motivation for 
the experimental work whose results will be tabulated 


ue 


it was conjectured by M. R 
those 


and discussed be low 


3. Experimental Data 


The algorithm described in the preceding section 


was carried out for the problem Ag 5 with 
wew ” 4 (h4erd ~ woe 
WA a Z24 4 ‘ 
Here ‘ or 4 41K ASS 
i 
oo ‘ ] ; lov 
2 ‘ 410 ‘ TT 
uM 270 al IT1IAe 4), 


and 


OoNeOY 
14279 
00024 
(M4578 
OO8O43 


OD4805 


Runs were made with fixed values of 8; ranging be- 
tween 6,=.1 and 8,=1.9. Since the matrix A and 


the vector 6 were obtained from an original ma 
B and vector ¢ as 


A=B*B, b—B*e, 


A is positive and symmetric. Hence, a conven 


metric 1s 


H=A 


In the metric (8) the gradient becomes &,—) 
Thus the other significant quantities defined in s 
tion 2 assume the form 


Table 1 lists runs of f(z,) and f(r)/f(x, fo 
various fixed values of 8;. In each case the imitia 
r=Z_9 was zero. The values of f(z,) are given her 
to the same number of places as in the original 
calculation, while the ratios f(x;)/f(z,_,) have bee: 
cut down from a six-place table. Table 2 gives th: 
values of a, corresponding to the runs listed in tab 
1. These numbers were originally computed to 10 
places. The results listed in tables 1 and 2 are based 
for the most part on single runs of the Card-Program- 
med Calculator. The main check used in the ma 
chine computations was the monotonicity of / 
However, the run 8,=.9 has been exactly duplicated 
on two separate occasions The results of the run 
8,=1 check closely with results obtained by A. | 
Forsythe in an entirely independent manner 
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The following 


We define the gradient as ¢ é 


Let 


gives the system of order 


i 


matrix whose le: 


Co 


found by 


st eigenvalue Is zero 


rresponding 
the 


YF 0 


MINI ZIng 


this 
original 


to 
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than or equal to one 
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By 


where re) 


Then 


y by the formula y 


al 
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view of 
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ues of B ranging between 
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A is 
components multiplied by 10 and e* 
O is me tually 
Runs were made 
7 and 


Starting 


riven 


Ly 


the 


algorithm was act 


by (5 


Sarre 


8] have examined in « 
convergence of various methods for a 


com 


letail 
plish 


algorithm 


ik that minimizes u 1s 
“optimum a’ procedure 
y Dy — wy ody 
pls x) 

IS a positive number less 
given y,, we determine 

y] at 
ally carried out with the 

h 

ene 
4 is given by (6) with all 
$335840 In 


a problem of the tvpe 


with various val- 
| For purposes 
point origin as 





in the computations recorded in section 3 was used 
and a record of 

y* Dy 

Yn+i 


was kept. In table 3° we present the results of the 
run § 9. The runs for other values of 8 showed 
substantially the same unstable behavior, while that 
for 8=1 (optimum) showed a very stable behavior 
and as a result converged quite slowly. The num- 
bers f(z,) and a,, appearing in table 3, were orig- 
inally computed to 10 places The ratios f(z ,)/f(z+~1) 
were computed from the original values of f(z,) and 
then cut down to the present size 


rarLte 3 


333840 

121482 0 
86700 

70323 

57533 

OOOUS 


3639 

7137 2575 
S111 7076 
SIS] 3951 
0550 5478 


S825 3968 
0057 5494 
9242 3479 
9420 5237 
“744 3046 


53564 
538640 
10788 
16002 
15704 
4020 5577 
9947 2761 
S853 6630 
9910 2482 
S607 9109 


11266 
11045 
36337 
30011 
31318 


0672 
8226 
S867 
3380 
9324 


2189 
7145 
1vOu 
9172 


2758 


30201 
24919 
22005 
7468 
6963 


= a 
Socmws 


B5U5 
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8378 
4060 
7540 


OOS! 
2493 
S282 
2201 
5481 


SUBS 
5541 
1642 
1206 
3171 


to bo BO Oo bo 
ee 


1Y12 
» 6970 
1766 
1380 
1008S 


2044 S334 
264 OOOR 
151 5733 
117 7708 
01 7806 


z 


to tt 


} 7 5 5 


1800 


11 1274 
be! 6933 


5. Conclusion 


The error function f(z) goes monotonically to zero 
in each of the gradient methods. Hence the number 
P,=100f(z,)/f(a) tells us what percentage of the 
distance from the starting point to zero remains to 
be covered at the nth step. Table 4 gives Py for 
various values of 8, We note 


* The author thanks R. M. Hayes for furnishing him with most of the data 
appearing in table 3 


TABLE 


8. 14 9. 05 9. 07 9. 10 


after 30 steps that, with the exception of the « 
8,=8=.1, all the gradient methods for values 
are converging faster than the optimum method 
while all the gradient methods for values 8 
are converging at a slightly slower rate than 
optimum method. For the eigenvalue method we 
have P..=.002. The increased rate of convergence 
for this case is offset on the Card-Programmed 
Calculator by the greater length of time needed for 
each step. However, on a high-speed machine this 
factor would be negligible. An explanation of the 
speedier convergence of the eigenvalue procedure 
lies in the fact that the transformation of the prob- 
lem to the homogeneous form has shrunk the ratio 
of the largest and the smallest nonzero eigenvalues 
This improvement of “condition” is something that 
one could not generally expect to occur [10] 

The ratios f(r, /f(z, compare the rate of con- 
vergence at each step with that of a geometri 
progression having the same ratio. A study of these 
numbers and of the corrections a, brings into sharp 
focus the contrast between the instability of the 
gradient methods employing 8<1 and the stability 
of those emploving 8>1. In the method of steepest 
descent it is just this stability that permits accelera- 
tion. However, this acceleration must be achieved 
through a modification of the computational routine 
On the other hand, the instability of the methods 
using 8<1 leads to occasional accelerations without 
the introduction of any changes whatsoever in the 
computing routine 

From the point of view of using the “almost 
optimum” gradient method on a large scale computer, 
its self-acceleration property has more than theoreti- 
cal interest. As is well known, the high-speed 
memory capacity of the computers now in existence 
is rather limited. Hence, the necessity of storing a 
special acceleration routine might prove to be a 
severe handicap indeed. 

It is worth while to compare the values of a, with 
the reciprocal eigenvalues of the matrices A and ID) 
For A these reciprocal eigenvalues have been found 
to be approximately 2.0, 3.9, 5.7, 12.1, 63, and 372 
while for D the finite reciprocal eigenvalues rang: 
between approximately 8.35 and 0.189... We wil 
pay particular attention to the points at which 
acceleration took place. One sees that preceding 
an acceleration there was a “smoothing run”’ during 
which the a,’s were in the range of the small recip- 
rocal eigenvalues. On the iteration immediatel) 
before an acceleration, a, was almost equal to the 
smallest reciprocal eigenvalue, while on the iteration 


‘These values were furnished by R. M. Hayes 


412 





during which acceleration took place a, was between 
the highest and next to highest reciprocal eigenvalues 
It was just this technique of choosing a, that the 
present author helped develop in previous experi- 
ments with a “fixed a’’ gradient method, which were 
onducted under the direction of M. R. Hestenes 
In this method the operator chooses the value of a 
just before each iteration, and by judicious choices 
he can successfully accelerate the method to a con- 
siderable extent. However, this requires too many 
judgments of the operator to be practical for a fast 
machine or an inexperienced operator. Hence, it is 


juite hopeful to note the existence of methods 
suitable for high-speed machines that can duplicate 


the fixed a@ acceleration procedures without any 
intervention by the operator once the process has 
started 
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An Alkaline Solution of Potassium Chromate as a 
Transmittancy Standard in the Ultraviolet 


Geraldine W. Haupt 


he reliability } photometric scale of spectro 
rlasses ar nsuitable) has led to the study 


chromate having th composition 0.0400 @ liter 


where 


ol an aqueous solution of .OTASS 
kK.CrO, in 0.05 NV KOH ased on extensive measurements ising photographic, photo 
electric, and visual spectro ometry, standard values of spectral transmittancy have 
been determined for 1.000 and 2.000 em of solution at 25° C from 220 to 500 my These 
values have been tabulated ong with derived values of absorbancy molar absorbaney 

er values have been compared graphically with those 

*hanges in temperature have 


ty, and 
ybservers he res in spectra 


transmittancy with cha 
iulso included a study of the effects due to the container 


termined 
er periods up to 8 years, and a comparison of results ol 


ther from k.( rQ), or tron K.CreO as one of the reage 


To make the solution most useful as a calibratio 
standard of spectral transmittancy in the ultraviolet 

\ means of testing the reliability of the photo- | information should be available regarding its per 
metric scale of spectrophotometers in the ultraviolet | manence, or stability. under specified conditions 
including the effects of the container bottles on th 
spectral transmittancy) and regarding any change in 
its spectral transmittancy with temperature over the 
range of room temperatures to be expected. Such 


l. Introduction 


region has been in demand for a number of vears 
With the introduction of commercial photoelectric 
spectrophotometers for the ultraviolet in 1940, the 
demand increased considerably 

Although glass filter standards of sper tral trans- data have accordingly been obtained 
This investigation was started in 1940 and discon 


mittance [{1, 2] ' have proved to be a very satisfac- 
tinued in 1943 because of the war. It was resumed 


tory and valuable means of testing the reliability of 
the photometric scale of spectrophotometers in the | in 1948, and preliminary data have been published 
visible region, they are in general unsatisfactory in | [8, 9, 10]. The present paper describes the work 
the ultraviolet. Not only do all but a few special | leading up to the determination of these preliminar 
glasses absorb completely below 290 my, but, in | values, based solely on the measurements made in 
addition, the ultraviolet spectral transmittances of | 1940-43, and describes the recent work on new solu- 
many glasses change on exposure to strong ultra- | tions with a new instrument, on the basis of which 
violet radiant flux [2, 3, 4]. In these respects solu- | the final values were derived as herein published 
tions in cells with quartz end plates are found to be 

superior to glasses for use in the ultraviolet In 2. Preparation of Solutions 
addition, solutions are reproducible and can be pre- 
pared in the laboratory from specifications with All the potassium chromate solutions studied were 
relatively small time and cost. Various solutions | prepared in the Bureau's Chemistry Division by 
have been studied or recommended for this pur- | W. Stanley Clabaugh. The solutions are of two 
pose. Preeminent among these is an aqueous solu- | general types (A) those prepared by dissolving th: 
tion of potassium chromate, K,CrO, (0.04 g/liter) in | salt in dilute solution of KOH, and (B) those prepared 
0.05 N KOH. This has alternate regions of high | by dissolving the salt in distilled water 
transmission and absorption in the ultraviolet, and The type (A), or alkaline, solutions were prepar 
this one solution in 1- and 2-cm thicknesses covers | as follows: 

the transmittancy scale from about 0.90 to about (A;): A solution of K,CrO, stock material, rea 
0.01 (0.046 to 2.0 in absorbancy). Several investi- | gent grade, 0.0400 g/liter, in 0.05 N KOH (solution 
gators have previously worked with this solution, | of potassium hydroxide prepared by dissolving 3.3 ; 
and from the information available in 1940 this ap- | of potassium hydroxide sticks (85% KOH) of reagent 
peared to be one of the best for the purpose. It was | ——"— ; othe ene 


! 


hese solutions have their limitations for use as cal 


accordingly selected for calibration.‘ the shortwave end of the visible region and in the ultraviolet region. | 

ena Sg AE visible, neither of them has sufficient absorption in the vi to be of muct 
for the purpose In the ultraviolet, the copper sulfate solution varies 1 
rapidly from high transparency at wavelengths above 3 my to high absor 
it wavelengths below 270 mg and thus has no general utili over the w 
region Che cobalt ammonium sulfate solution used in the Davis-Gibso 
7} has too high tr ney between 230 and 400 ma t of much value 
spectrophotometric standa Solutions prepared at 10 times this concentrat 
were found (unpublished « to be unstable with time throughout the 
violet and were thus considered undesirable for use as spectrophotor 

ilibration standards 





quality in sufficient distilled water to make 1 liter 

A A solution of K.CrO, of the same concen- 
tration and alkalinity as (A,) but prepared from 
g of KoCr,O,/, which, when converted gave 
0.0400 ¢ of K,CrO,/liter This follows from the 
reaction K,Cr.O,+2KOH =2K,CrO,+ H,O 

The purpose of dissolving the potassium chromate 
n dilute solutions of KOH was to prevent any di 
hromate from forming. In A, a potassium chro- 
nate solution was prepared by means of a different 
This solution was 


O O30 


material, potassium dichromate 


epared to determine whether the spectrophoto- 


tre data from A, and A, would be identical 
Since thes found to be identical, there is a 
tinct advantage in using potassium dichromate as 
of the reagents rather than potassium chromate, 
cause in stock material potassium dichromate exists 
a pure! A special sample of potassium 
chromate that had been recrystallized four times was 
the preparation of one of the B, or neutral, 
tions with the view to studying the effect of the 

tv of the potassium chromate 
The solutions prepared in 1940, 1942, and again in 
10449 were stored in liter bottles of the ordinary 
storeroom glass type with glass stoppers. For fur- 
ther study of the effect of storage, half of the 1940 
alkaline solution prepared from K,CrO, and half of 
the 1940 neutral solutions were stored in similar 
bottles lined with ceresin The tips of the rubber 
stoppers used were also covered with ceresin. Also 
half of the alkaline solutions pre pared by methods A 
A, in 1949 were stored in alkali-resistant glass 


Corning glass 728 


were 


state 


None of the solutions was exposed to light except 
when actually in use to fill the cells for their measure- 
ment No solution was ever poured out of its bottle 
In filling the cells a pipette was used to draw off the 
Where flaking occurred, the sediment thus 
remained on the bottom of the bottle. Each cell 
containing solution measured relative to a 
similar cell containing distilled water 

The alkaline solutions are those from which the 
standard set of spectral transmittancies given in 
table 1 were obtained. Although both alkaline and 
neutral solutions were prepared originally for the 
purpose of establishing standard spectral transmit- 
data, the neutral solutions were finally not 
used for that purpose. However, as the neutral 
solutions have been studied over a period of 8 years, 
a summary of the results on these is included in 
section 8, 


3. Methods of Measurement 


The solutions were measured on five instruments 

three different methods. One instrument oper- 
tes by a photographic method, primarily for the 
traviolet region, but extending also into the violet 
d blue; one by a visual method; and three by 
\0toelectric, one of these being for the visible range 


solution 


was 


tancy 


obtained from the 
e prepared, the 


National 


chemicals 


3.1. Photographic Method 


The photographie data were obtained by means of 
the Hilger sector photometer with the Fuess quartz 
spectrograph [11, 12, 13] 

The equivalent slit widths for an 
which throughout all the 


0.2-mm slit 


was used exposures, are 


On each plate exposures were made of at least two 
aluminum spark spectra for wavelength calibration 
\ calibration curve was used to translate 
scalar values, which were read at the density match- 
points, into wavelengths Both 1.000- and 4.000-em 
cells were used with the Hilger instrument 


purposes 


3.2. Visual Method 


The Kénig—Martens spectrophotometer [13, 14] is 
a visual instrument. The slit widths used on this 
instrument are varied, depending upon the wave 
length and the spectral characteristics of the sample 
being measured. With an incandescent source the 
equivalent slit widths in millimicrons used for two 
wavelengths and for two settings of the slits 
riven below 


are 


Wavelength 


The mercury arc, also available for use on this 
instrument, enables transmittancy measurements to 
be made that are not subject to possible wavelength 
or slit-width errors 

Two cell lengths were used with 
namely, 4.000 and 10.000 em 


this instrument, 


3.3. Photoelectric Method 


Double Monochromator and Auxiliary 


Equipment [15 


a. Brackett Quartz 


The source of illumination used was a special incan- 
descent lamp, with a quartz window, which operated 
on a 120-volt, 45-ampere battery. The monochro- 
mator has three slit systems. The exit slit next to 
the photoelectric tube was kept at a constant mechan- 
ical width for nearly all the measurements, and the 
other two were varied concurrently Any stray 
energy that passes the second slit is effectively elim- 
nated from reaching the exit slit by the use of the 
second spectrometer provided on this instrument 
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The instrument operates with a single wavelength 
Spectral bands transmitted by the slits 
were comparable in width (my) to those obtained 
with the photographic method. A photoelectric cell 
of the electron-emission type was used as detector 
and values of spectral transmittancy were obtained 
by means of the ratio-of-deflections method with a 
high-sensitivity galvanometer 

The seale for reading or setting wavelengths on this 
monochromator is a linear one (not direct reading 
in millimicrons) and a wavelength calibration curve 
of the instrument was already available. This 
curve was based upon scalar settings previously deter- 
mined for mercury lines down to 237.8 my. As the 
incandescent lamp was the only source of irradiation 
for the spectrophotometric work, the wavelength 
calibration was frequently checked at one wavelength 
557 mu) by means of two particular standard glasses 
whose spectral transmittances are known to be the 
same at that wavelength. The linear correction 
determined at this wavelength was then applied to 
the whole scale 

The lower end of the wavelength calibration curve 
was extrapolated below 237.8 my for these measure- 
ments, and the data in this region with the Brackett 
instrument are therefore of lesser validity 


control 


b. General Electric Recording Spectrophotometer 


The General Electric recording spectrophotometer 
has been described in a series of three papers [16, 
17, 18]. The data obtained with this instrument 
are for the visible range only and extend from 400 to 
750 mu. Checks on the wavelength scale were made 
in the usual manner by using a calibrated didymium 
glass. A nominal slit width equivalent to approxi- 
mately 4 my was used throughout the recordings 
Cells of 1.000-cm thickness were used for the meas- 
urements 


c. Beckman Quartz Spectrophotometer (Mode! DU) 


The Beckman spectrophotometer 
scribed in a paper by Cary and Beckman [19] 
Transmission measurements with a modification of 
this instrument, which has a constant-temperature 
enclosure designed and constructed at the National 
Bureau of Standards, have been described in detail 
in a paper by Gibson and Balcom [20]. All measure- 
ments on the Beckman spectrophotometer were made 
with a 1.000-cm cell. The cell holder is kept constant 
at 25° C or other desired temperature by means of 
water pumped through from a constant-temperature 
bath. 

For wavelength checks the Hg 
mu) was used. When necessary, 
length scale was adjusted by the 


has been de- 


green line 546.1 
the whole wave- 
slight turning of 


* Acknowledgment is here made to W. F. Roeser and M. 8. Van Dusen for 
naking the instrument and its calibration curve available 

’ In 1942 a few measurements on a Beckman quartz spectrophotometer in the 
Bureau's Chemistry Division were made on an alkaline solution prepared in 
140. No ] ‘ ble for these me rement 


temperature contre Was aVaila sure 























the internal screw that rotated the collimat 
mirror. After such check or adjustment, the pre: 
wavelength correction was then applied at e; 


wavelength setting in accordance with the cu 
illustrated in figure 6 of reference [20] 
For transmittancy measurements three sour 


are available, namely, the incandescent-filam: 
lamp, the hydrogen are, and the mercury arc, all as 
supplied by the maker. The incandescent a 
hydrogen sources were used over the appropri: 
wavelength ranges, overlapping from 320 to 380 n 
The slit widths with these sources never exceed 
2 my and were usually less than 1 mp. The mereu 
arc was also used for transmittancy measurements 
at certain wavelengths as noted below 

A clear quartz lens similar to that supplied with 
the instrument was used to effect an approximate 
collimation of the beam. Over the wavelength 
range from 330 to 420 my with the incandescent 
source and from 360 to 380 my with the hydrogen 
source, a purple (Corning filter No. 9863) lens was 
used to reduce stray light as well as to collimate th 
beam. 


4. Standard Spectrophotometric Data 


In examining the data obtained during the first 
two vearsof the investigation, the type A, or alkalin 
solutions were found to give more reproducible re- 
sults (i. e., less scattering of data on repeated meas- 
urements) than the type B, or neutral, solutions 
Furthermore, the alkaline solutions were found to 
give a better range of transmittancies in the extreme 
ultraviolet than the neutral solutions. Since the alka- 
line solutions were thus found to fulfill more nearly 
the requirements of a standard of spectral trans 
mittanecy in the ultraviolet, all later work and all 
the detailed data to be presented refer to the alkalins 
solution of potassium chromate (0.04 g/liter in 0.05 
N KOH). 

As already noted, the work on the 1940 and 1942 
solutions was interrupted in 1943. Sufficient data 
were obtained, however, that a tentative set of 
spectral transmittancy values could be evaluated 
These values are given * in column 2 of table | 
designated as 7’. They were read from a curv 
of which figure 1 is a tracing. The individual deter 
minations by the five instruments used are plotted 
and may be identified from the legend. Thes: 
tentative values thus represent a graphical average 
of all these earlier data. Measurements were mac 
on the 1940 solutions a few days after their prepara- 
tion and then again at the end of 2, 4, 7, 8, and 10 
months, 2 years 4 months, 2 years 5 months, and 
2 years 7 months. The 1942 solutions were meas- 
ured a few days after their preparation and again 
in 1, 2, and 6 months. 


See 


* These tentative values have‘also been published in [8; 9, p. 43-44; 10, p. 2 
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In (April 1949 four sets of data were taken on the 
20), on each of 
freshly prepared alkaline \, and \, 


tained at 25° C both in glass and 


Jeckman spectrophotometet one 
solutions main- 


stored alkali- 


esistant ware 
Che four sets of data were averaged Values ob- 
tained with the mercury source are listed in column 


the hydrogen 
in column 4 
were 


of table 1; values obtained with 


are listed 
in column 5 


1 ol incandescent sources 
The final values of 7 

ved according to the following procedures 
At the a straight aver 


if three sets of data 1949 Beckman 


show ni 


mercury wavelengths, 


namely ] 


a with mercury line source (column 3); (2) 1949 
Beckman data with hvdroge n and/or incandescent 
nuous-spectrum source column a: & Tenta- 


column 2 


with the 
of two 
with 


ies based on 12 data, 7°, 
\t other wavelengths above 220 mu 
295 to 315) a straight average 
namely 1) 1949 Beckman data 
ent column 4 


data, 7 


1940 


on ol 
data 
souree 


1940-42 


een and or incandes 


Pentative values based on 


) 


olumn 
At 295 and 315 may slight further 

stments in 7 exceeding 0.003) were made 
final values into better agreement with 
wavelengths, 


S00. 305. 310 


hot 
to bring the 


he averages obtained at the mercury 
ture (a 


and S of tabl 


y proc 
| show the deviations 
he three sets of data from the final values It is 
nt that there mav be considerable uncertainty 


he third decimal of column 5 


Columns 6, 7, 


Values of spectral absorbancy, Ay, for a 1.000-cm 
kness of solution are given in column 9. These 
is of A, were derived from the relation 

| logy, 1 To cover the photometric scale be- 
/ 0.10, values of transmittanev for a 2.000-em 


were computed and are given in column 10 
Column 11 a the molar-absorbane\ 
21] ay 1./be where ¢ is the concentration 
| ressed in moles pet liter and h Is the cell thick- 
Column 12 gives the values of 


es values of 


rh x 


ss in centimeters 


5. Comparison of Data with Those of Other 
Observers 


Tabular data on the spectral transmittancy of 
potassium chromate in 0.05 N potassium hydroxide 
have been published by von Halban and Siedentopf 
22}, Réssler [23], Baly, Morton, and Riding [24] 
ind Hogness, Zscheile, and Sidwell 25) 

Hogness, Zscheile, and Sidwell used a photoelectric 
method for their determinations. The potassium 
hromate was recrystallized four and five times. 
They show the absorption spectrum for this solution 

a figure [25]. The dots through which the curve 

drawn in their figure are taken from several series 

determinations at 2.5-mu intervals from 220 to 
110 mu. Values of molar extinction for the maxima 
nd minima and for 366 and 265 mu were given in 
bular form. The experimental conditions under 
hich these tabulated data were obtained are givep 


table 2 


Figure 2 the data with those 
already published in tabular form 
for this comparison in accord with the practice used 


in the International Critical Tables [26]. The meth- 


compares present 


Log (ly, 1S used 


ods used in the various studies are shown in the 

caption 

PaBLe 2. Experimental cond ’ hich certain data 
of Hloqness ‘scheile. and Nidwe r ere obtained 


seen that in general the 


Krom figure 2 it may be 
agreement In logyd@y from 220 to 450 my found by 
the five investigators shown is quite close, with the 


exception of a few points. Below 300 my, two values 


of von Halban photoelectric method) at 254 and 
°265 mu are at variance with the rest Above 300 
my, the value of Baly, Morton, and Riding (photo 
graphic method) at 313.2 my, seems erroneous 
Most of the 16 wavelength values of Réssler from 
381.5 to 439.0 mu will be seen to be lower in logyay 
than those of the others shown. In his paper 


Réssler has put these wavelengths in parentheses; 
apparently indicated the 
relatively low dispersion of the prism in his instru 


these were so hee ause of 
ment 

In the present investigation only a single concen 
tration of potassium was No 
attempt was made to carry out a study of conformity 


chromate used 


to Beer’s law, because the one solution in 1- and 
2-cm thicknesses covered the photometric range of 
interest However, it was known that Hogness, 
Zscheile, and Sidwell had made a run on exactly 
the same concentration used by von Halban and 
found good agreement with their curve [25] obtained 
from data on optimum concentration. From these 
data it was concluded that Beer's law held Results 
have also been published 27) of Beer's law tests on 
potassium chromate in 0.05 N KOH at 372 my 


Various dilutions were used, and the results are given 
in both tabular and graphical forms. The molecular 
concentrations used varied from 1.61% 107° to 51.5% 
10 This corresponds to vamations in concentra- 
tion from 0.00313 to 0.1000 g/liter. Determination of 
this compliance with Beer s law was made asa means 
of checking the photometri reproduc ibilitv of a new 
spectrophotometer 


6. Effect of Temperature on Spectral 
Transmittancy 


Because it was desirable to know whether or not 


the potassium chromate solution for use as a calibra- 


tion standard changed in spectral transmittanecy 
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Fiaure 2 


Molar solution in 0.05 


Values of log oO for potassium chromate. 
N KOH. b=1 cm 


au 1,/be = —logy T./he. 


s. From NBS final data, column 12 of table 1 


b. Hogness, Zscheile, 
ec. Réssler; photographic 


d. Baly, Morton, and Riding-Hilger rotating sector 
and Riding Judd-Lewis sector 


Baly, Morton 


von Halban; photoelectric 


with change in temperature, several preliminary 
tests were made in 1948 on the 1942 solutions by 
means of the Beckman spectrophotometer [20]. The 
resulting data indicated that the largest effects due 
to changes in temperature occur below 235 mu 

A more detailed and complete study of the change 
in spectral transmittancy with temperature was then 
made for the 1949 alkaline solutions of both types 
stored in alkali-resistant bottles. The wavelength 
range from 215 to 500 my was studied for approx- 
imately a 25-deg change in temperature. The 
spectral-transmittancy data were determined for one 
solution over a definite wavelength range at 25°, 50°, 
and again at 25° C in a single day. Each change in 
spectral transmittancy for an approximately 25- 
deg change in temperature was computed separately 


for each day’s set of measurements. Three sets 


Sidwell; photoelectric 


photographic 
photogr aphic 


were averaged and are plotted in large circles in 
figure 3. 

Included in the figure are also points obtained 
from the measurements made on the 1942 alkaline 
solutions at various times at each of the tempera- 
tures, 25° and 50° C. These data corroborate the 
data obtained with the 1949 solutions. It may be of 
interest to point out that the changes in transmit- 
tancy with increase in temperature shown in figure 3 
are such as to indicate a shift of the absorption 
bands toward longer wavelengths. 

The curve shown in figure 3 was then drawn, 
from which values of A7, for a change of 25 deg 
were read. Values of AT, for a change of 10 deg 
were computed from these curve readings by linear 
interpolation, and the resulting data are given in‘ 
table 3. Users of the standard solution should 


420 
































. 
> 
. 
.° 
. 
> 
. 
° . 
‘ . 
o 
> 
. 
° Pa 
. 
id %. . 
2 el 
2 ‘e ot + tame* _ 
° _— 
) 
j 
| . 
| . 
| 
| > 
} . . a 
| . A) . @) 
9 p 
} . 
| . * 
Q . A } . *-e . fe 
. * 
. oO o- , . — “oO 
| . . 
| ee . . 
\ ‘ . . oO . 
| - * Ow »s ° " 
oO 
of oO 
} o . ~ . 4 
| . . . 
i o J 
| By 
. } 
5 
| } 
| > 
| } 
4 
. } 
A r be r 
WAVELEN [4 
Ficure 3 Change in spectral transmittancy for a change in temperature of 25° C 
In establishing this curve more weight has been given to the data obtained from the 1949 solutions The values of AT, given in table 3 for a 10° change in temperature 
were derived from this curve by linear interpolation 
>A Ti( = T: Ts) for 1949 solutions he measurements were all made on the same day 


@.AT, for 1942 solutions. (The measurements were made in 1948, mostly on different days 


inake use of these corrections proportionately 
whenever the temperature of the solution deviates 
importantly from 25° C. Values of AA, in table 3 
are derived from the values of A7, by the relation 


AT, 


s (ave 


AA, logo ¢ 


7. Permanence 


Measurements have been made from time to time 
these alkaline solutions, and the study has in- 
ided effects due to the storage of some of the 
lutions over periods up to 8 years. The 1940 
ukaline solution made from K,CrO, in KOH was 
tored both in an ordinary glass bottle and in a 
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similar bottle lined with ceresin. The 1940 alkaline 
solution made from K,Cr,O,; and the 1942 alkaline 
solutions made both from chromate and dichromate 
were stored only in ordinary glass bottles. 

Above 260 my no definite change in spectral trans- 
mittancy was found for the solutions stored in glass 
Below 260 mu the 1940 alkaline solutions showed 
definite changes at the end of 8 vears of storage 
in glass bottles, the solution made from the dichrom- 
ate exhibiting the smaller change. The fact that 
potassium dichromate in stock material exists in a 
purer state may account for this solution being more 
stable. However, the 1942 alkaline solutions of 
both types likewise stored in glass bottles showed 
little if any change in spectral transmittancy over a 
period of 6 years. 











TABLE 3 Value \7 7 Tx, or Ty ] 1 10°C 
change ir ’ wresponding to wa noth 
7 
] 
1 
‘ Cin} 
aot * 
' . , 
7 a aM 
¥ 7 x mM 
, ‘ i 
, wn , 0s 
, wi a , - 
For the 1940 alkaline solution made from K.CrO, 


in KOH and stored in a ceresin-lined bottle, measure- 


ments from time to time showed that changes in 
absorption were occurring. Data obtained at the 
end of an *&-Vveal period showed that complete 
absorption had occurred below 370 mu. Obviously, 


therefore, an ordinary glass bottle is better 
storage of a solution of potassium chromate in KOH 
than a ceresin-lined bottle. No study of the per- 
manence of the alkaline solutions stored 
resistant bottles has been made 

It must concluded from that 
since the data below 260 mu show possibility of a 
change in spectral transmittancy after a 6-month 
period solutions ho older than 6 months should be 
used for the photometric calibration of a 
photometer in this region of the spectrum 
other hand, for regions above 260 mu, it seems safe 
use alkaline solutions of either stored in 
ordinary glass bottles up to a period of 5 vears 


in alkali- 


also be this study 


spectro- 
On the 


to type 


8. Neutral Aqueous Solutions of Potassium 
Chromate 


Considerable data were obtained on the B solutions 
(see section 2) and a summary is presented here 

The spectral-transmittancy curve of neutral solu- 
tions of potassium chromate was found by Hantzsch 
[28] to differ somewhat from the alkaline solution. 
and the results reported here corroborate this finding 
From 400 to 500 mu the curve for the neutral solution 
is less steep than that of the alkaline solution with a 


crossing point at about 440 mu, above which the 


for 








transmittancy of the alkaline solution is slig! 
higher. Between 300 and 330 mu the neu 
solution is slightly lower in transmittancy than 


alkaline solution and slightly higher both below 
above this range up to 440 mu Below 230 mau 1 
difference is large, the transmittancy at 215 
being 0.40 for the neutral solution and 0.07 for 
alkaline The value (0.07) for this wavelength n 
be found in table 1, column 4. More weight sho 
be given the 1949 Beckman 225 
because the measurements 
closely controlled conditions of temperature 
while the solutions were still fresh, unport: 
controls for this particular region (sections 6 and 
However, this particular value may be somewhat 
the 


data below 


were obtained un 


I 


two 


error due to another cause, namely 
presence of impurities in the KOH 
Over a period of years the neutral solutions we 
found to be permanent than the alkali 
solutions. Significant 


poss! 





In a giass bottle were confined to the spectral regu 
below 230 mu: and for storage in a ceresin-lu 
bottle changes were found only below 260 mu. A 
specially pure sample (reerystallized four times 


changed still less. The neutral solutions 
cause flaking in the bottles as the alkaline solutior 
do 

From the standpoint of 


venience in storage the neutral solutions are slight 


permanence and cot 


superior to the alkaline solutions, and had they beet 


found equally reproducible the neutral solution 
might be preferable spectral 


standards instead of the alkaline solutions 


as 


9. Summary 


Spectrophotometric studies of neutral and alkali 
aqueous solutions of potassium chromate extendil 
over a period of vears have shown that the alkali 
solution is the better suited for spectra 
transmittancy standard in the ultraviolet. TI 
spectral transmittancy of this solution has been « 


use as a 


termined bv five spectrophotometers used according 


to three different methods, and the final values : 
given in table 1 from 220 to 500 mug. The final 
sults are also expressed there in terms of absorban 


changes from S8-vear storag: 


did hot 


transmittancy 













| 


A log.7,, for a l-em cell; molar absorbances 
index, vg Ale he . and logy vy Changes of trans 
mittancy with temperature were evaluated, al 


correction data for deviations from the standard ten 
perature of 25° C are given in table 3 

When stored in ordinary glass bottles, the solutior 
is reasonably stable for 5 or 6 vears. For work of t! 
highest precision, however, avoid unsuspect 
errors from any cause, it is recommended that so 
utions not over 6 months old prepared from chemi 
of the highest purity be used as transmittancy sta! 
ards. 


to 
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Two Applications of Group Characters to the Solution 
of Boundary-Value Problems 


E. Stiefel ’ 


It is shown that the numerical work involved in solving a boundary-value or eigen- 


value problem by finite difference methods in a domain with many symmetries can often 
be reduced by applying the theory of group characters to the group of symmetries of the 


domain 


In the second part the author considers the problem of solving Au=0 in a cube when 


the prescribed boundary 


values are invariant under the group of rotations of the cube 


In the series representation of the solution in terms of harmonic polynomials only a subset 
of these polynomials actually occurs, and the theory of group characters facilitates con- 


siderably the determination of this subset 

It is very well known that representation of groups 
may be a useful tool in discussing various eigenvalue 
problems in quantum mechanics. This is due to the 
fact that boundary-value problems of Schroedinger's 
equation are to be solved mostly in the whole space 
or inside a sphere as fundamental region and therefore 
the high symmetry of those domains can be used. 

It is the purpose of this paper to show that group 
characters sometimes are appropriate also to sim- 
plify the numerical computation of boundary-value 
or eigenvalue problems in more general-shaped do- 
mains, in particular if the problem has been trans- 
lated into the language of difference calculus It 
will be sufficient to explain this discussing simple 
examples 

Let, for instance, A be any plane region and D a 
differential operator invariant under Euclidean mo- 
tions or reflections of the plane. More exactly, if 
J(z,y) 1s a function in our plane and .V/ a Euclidean 
rotation, translation, or reflection, then 

D( Mf) = M(Df), (1) 
which is to say that D must commute with any 
motion M. The Laplace operator A and the operator 
AA of elasticity problems have this property. 

In order to solve any given boundary-value prob- 
lem with respect to the region A and the operator 
D, we may use the difference technique, introducing 
a square lattice L in A and replacing / by a lattice 
function (defined only in the lattice points) and D 
by a difference operator, which approximates the 
given differential operator. In the case of a A-prob- 
lem, the corresponding difference operator may be 
given by figure 1. Relation (1) will be true again in 
the difference case, provided that by ./ is understood 
any motion or reflection transforming the lattice 
into itself. 

Suppose now that the given region A is invariant 
under a group G of transformations of this type. If 
A has no symmetry at all, we may use Schwarz’s 
alternating process dividing A into overlapping sub- 
regions, each of them having symmetry properties. 
as niece Hoan 

This work was performed on a National Bureau of Standards contract with 
the University of California, Los Angeles, Calif., and was sponsored (in part) 
by the Office of Naval Research 


? University of California, Los Angeles, 
nische Hochschule, Zurich, Switzerland. 


Calif., and Eidgenossische Tech- 


FiaurRe 1 


It is our goal to show how the group G can be used ir 
order to simplify the solution of the given boundary- 
value problem in the region A 

Take, for instance, the square region A of figure 2 
The group ( is in this case the symmetry group of 
the square of order 8 having the following similarity 
classes: 


Number 


Class Elements of 
elements 
[1] Identity l : 
[2] Rotations +90 2 - 
[3] Rotation 180 l 
[4] Reflection in horizontal or 
vertical line 2 
{5} teflections in diagonal lines 2 
Total 8 


The given boundary-value problem is equivalent to 
a system of nine linear equations for the values of the 
wanted function in the nine interior lattice points 
Any lattice function f may be visualized as a vector 
in a nine-dimensional vector space S. Therefore, 
applying the operations of the table to the lattice 
functions, we get a linear representation FR of the 
basic group G in S. 
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FIGURE 2 


It is not too hard to show that its character x is 
given by the following rule: The value x(.MVJ) of x for 
an element .V of G is equal to the number of lattice 
points fixed under the operation \/. (Use as coordi- 

system in S the lattice functions not equal to 0 

in one single lattice point; they are permuted by 
\f.) We want to decompose the representation 2? 
nto irreducible components 


R=e,R,+eoR + ess, (3 


Where R,, Rs, ‘ R are the irreducible repre- 
sentations of G. (There are exactly five because @ 
has five similarity classes According to the ortho- 
gonality theorems of representation theory this has 
to be done working on the table of characters of @ 
as follows: 


(lass 1 2 3 { 15] 
No. ofelements | 2 ] 2 2 
x l 3 . 
X, l ! | l ] 0 
X l l l ] l I 
X, l l l l l l 
x 9 0 2 0 0 2 
x iT) l l 3 3 


In the first lme we have listed the similarity classe 
with their numbers of elements and in the following 
lines the values of the characters of the irreducible 
representations in those classes. They are recorded 
in many textbooks of representation theory.’ The 
last line is the character x of our representation P, 
computed according to the rule previously mentioned. 
The first number in each line is the degree of the 
representation (dimension of the vector space trans- 
formed by the representation). According to a 
general formula of character theory the coefficient 
¢, in (3) is 


1 
SS heed 


>) x(x MW), (5) 
Vl 


\. Speiser, Die Theorie der Gruppen von endlicher Ordnung, 3d ed., p. 179- 
Springer, Berlin; and Dover Put 


with the group element \/ running through the whole 
group G. By this we have for instance from our table 


The multiplicities ¢; are listed in the column to the 
right of the table. Finally, we find 


R 3R;, t R, t R, T 2R 6) 
Therefore, the vector space S splits into subspaces 


S 2 ks: & & & & 7) 
each of them being invariant under 7? and being 
transformed by the irreducible representation given 
by its subscript 

Take now a set of subspaces hav ing the same sub- 
script (i. e., being transformed by the same repre- 
sentation) and introduce base vectors in each space 
of the set. For instance, S; and S,’, being both 
two-dimensional, may be spanned by the base vectors 


Re Ge ¢ 


Generally speaking, each set of equally transformed 
subspaces yields a rectangular arrangement of base 
vectors by listing the base vectors of each subspace 
in a line. The length of the rectangle is the degree 
of the irreducible representation under consideration 
and its height the multiplicity of this representation 
as a component of ?. It is, of course, possible to 
choose the base vectors in such a way that each line 
of the rectangle is transformed exactly in the same 
Ww ay 

It follows then from a well-known generalization 
of Schur’s lemma that each linear operator in S com- 
mutable with the transformations of /? transforms 
the vectors of each column in our rectangle among 
themselves and in each column in the same way. 
Taking into account the relation (1), we see that the 
difference operator of the given boundary-value 
problem has this property. From this follows imme- 
diately that the original system of nine linear equa- 
tions with nine unknowns splits into partial systems 
The rectangle (8), for instance, vields two systems of 
two equations each, having the same coefficient 
matrices. From the first and the ¢ column of the 
table (4) we get therefore the following final result 
The system of nine equations splits into one system of 
three simultaneous equations, one single equation, 
one single equation, and two systems of two simul- 
taneous equations. 

If we have to deal with an eigenvalue problem in- 
stead of a boundary-value problem, the analogous 
statement is, that the wanted characteristic poly- 
nomial is the product of a factor of third degree, two 
linear factors and the square of a quadratic factor. 
In order to carry out the decomposition of the system 
of nine equations explicitly, it is, of course, necessary 


‘ 


‘B. L. van der Waerden, Die gruppentheoretische Methode in der Quanten- 
mechanik, sect. 13, p. 47-50 (Springer, Berlin, 1932 


425 














to establish the lattice functions, which build the 
base vectors of the different rectangles s The 
three bast lattice functions corresponding to the 
first irreducible representation /?; are simply three 
linea! independent functions invariant under all 
motions ol the square As another exampl we give 
the values of the lattice functions ‘ in the nine 
tte pomts 
) i) ) ‘) | {) 
| () ] é () () () 
Oo Oo 0 0 0 
y 
] 0) | | () | 
oO Oo O é ) ft 0 
] () | | () ] 
Taking as J) the operator of figure 1, we get 
Dh he é 
10 
De J be 


are transformed in the same wa\ 
D with 
equation 


The functions ¢ ‘ 
The eigenvalue problem of the 
vanishing boundary values gives the 


oper ntor 


Df+xrf=0 1] 


Let us put f= ae he Then (10) vields the eigen- 
value problem 


A\—4)a+2b=—0 
a h—4)46—0 
And from this follow the eigenvalues 
\= 44/2 13 


of the operator D, each of them having at least the 
multiplicity 2 

It is not too hard to establish for a lattice in a 
riven region A the rectangles of base vectors when 
we know the characters and the irreducible repre- 
sentations of the symmetry group of A. It 
venient to choose all those vectors (their total number 
is equal to the number of lattice points) orthogonal 
to each othe This is possible because the re- 
presentation # is orthogonal, and therefore the ir- 
reducible subspaces of S are orthogonal 

At the Institute of Applied Mathematics in Zurich 
(Switzerland) we used the methods of this paper in 
order to solve the boundary-value problem of 
Airy’s elasticity function inside a As subre- 
Schwarz's alternating 
squares of 16 lattice points 

In his book about eigenvalue problems, L. Collatz°® 
introduces the notion of “eigenvalues of a graph” 
and discusses an example, where the graph is built 
by the joining midpoints of 
edges of a cube. This problem may also be treated 


' Collet: } wertaufeaben mit t ’ ‘ 
Leipe 1Wie 


is con- 


dam 


rions ih process, we chose 


segments two 
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by characters, using the group of all rotations 
the cube 


In this section another application of chara 


theory is discussed, related to the methods of 
Bergman and \l Picone for the solution 
boundary-value problems. For our purposes th 


the 
a complete set of 


methods may be characterized by follow 
procedure. At first partic 
solutions of the given partial differential equat 

is constructed and then the solution of the bounda 
value problem is expanded into a series having as 
terms those particular integrals. Let us take 
instance, Laplace’s equation Au=0 im three-dimy 
sional space. The particular solutions may be 
this case the harmonic polynomials of degree n 


Y (2.9, 2 preg P“™ (cos 6), n=0,1,2 

14) 
m running from n) to (+n On the right side 
r, 6, @ are polar coordinates, /’,” is a generalized 
Legendre polynomial. Suppose now that the fun- 


damental region A of the boundary-value problem 
is a cube having its center at the origin and suppose 
furthermore that the given boundary values ar 
invariant under the group @ of the 24 rotations of 
the cube. It is obvious that under those cireum- 
stances only harmonic polynomials appear in the 
Bergman-Picone expansion, which have the samy 
symmetry properties. In other words, we have 
the problem to establish all harmonic polynomials 
invariant under the group ( of rotations of the cubx 

The theory of characters gives general methods 
in order to solve problems of this type. @ has th 
following similarity classes 


Number 


(‘lass rypical element of ele- 
ments 
1} Identity l 
2 Rotation O0 6 
3 Rotation 180° (z-axis 3 
} totation ISO” (axis: 2 Or y 6 
) Rotation + 120° (axis = diagonal S 
rotal 24 


By any rotation of the cube the (2n+1) harmo 
polynomials of a given degree n are transformed 
linearly among themselves and vield therefore 4 
representation FP, of degree (2n+1) of G Let M 
be any rotation of G and x,(.\/) the character of / 
The number of linear independent harmonic poly- 
nomials of degree n invariant under @ is equal to 
the multiplicity ¢c, of the unit representation in th 
decomposition of FR, into irreducible componen! 
(The unit representation of a group maps eve! 


This question has been posed and answered by G. Pols and B. M 
Sur les svmé@tries des fonctions sphériques de Laplace, pt. rend, 228 0 
1949 Sur les fonctions sphériques ae Laplace vmétrie cristallogra 


lonnée, Compt Te 1049 

















roup element into the one-row matrix= 1 This 


multiplicity is given by the character formula 


e,.=—~ DS *x,( MM 16 


inning through the group. It remains to com- 


the character x V/ In order to do. this 
mav restrict ourselves to rotations .V around 
he c-axis. Let gw be the rotation angle. Taking 


ccount (14), it follows that the (2n+-1) har- 
polynomials of degree » are transformed by 
the following wavy 
y * r y 17 
matrix of the transformation is pure diagonal 
the trace 


Krom (16 and the table 15) follows now the result 


ae sin (27 1)45 
( 
24¢ sin 45 


sin (2n 1)60 
9 sin (2n+1)90 8 = 19 
sin 60 \ 


The final result may be expressed in the following 
terms. Write » in the form of a multiple of 12 and 
the corresponding remainads I 


12-k+-) O<; ) 2) 


Then 


where ¢(r) is given by the table 


The poly nomials themselves can be built by « hoosing 
any harmonic polynomial of degree n, applying the 
24 rotations of the cube and adding up the 24 result- 
ing polynomials 

A more general probl m would be to establish the 
harmonic polynomials being transformed by @ 
according to a given irreducible representation ol G 
This can be used to simplify the solution of a bound- 
ary-value problem in which the given boundary 
values have not the symmetry of the cube. The 
basic idea is, of course, the same as outlined previously 
in this paper 


Los ANGELES, January 29, 1952 
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Heterogeneous Equilibria at the Glass Electrode-Solutio: 
Interface 


Donald Hubbard and Richard G. Goldman 


Heterogeneous equilibria at tl 
with the pH response of glasses 
the glass surface by acid aqueou 
gratable negatively charged ior 

system in accordance with tl 
of products.”’ The unever 
conditions was readily 
ia series of Na,O-SiO, glasses 
the various experimentally 
of swelling by high co 


} 


i 


SiQ>, an 


repressior 
and hydrofluoric-acid solutions 
metry potentials However, the 
preferential voltage response of the 
partial exclusion of other ions, 


l. Introduction 


upon immersion in most acid 
aqueous solutions exhibit a preferential leaching 
of alkali and other cations from the more durable 
framework. This remaining silica-rich layer acting 
as negatively charged nonmigratable must 
cause an distribution of all 


rlasse 
lasses 


Silicate 


Ons 


necessarily uneven 


migratable ions in the system in accord with the 


distribution law and Donnan membrane equilib- 
rium [1, 2].'. Thus, a glass surface, which is generally 
looked upon as chemically inert and static, turns 
out to be transient for every solution and seething 
with activity. However, the framework of the 
silica-rich layer for many glasses is frequently 
sufficiently resistant to the attack of a large number 
of aqueous solutions to establish conditions that 
are adequately stable to be studied and discussed 
276] The present 

ascertain what 


in terms of equilibria [3, p 
investigation was undertaken to 
experimental evidence could be obtained demon- 
strating the nature of the heterogeneous equilibria 
at glass-aqueous solution interfaces and to determine 
what part these equilibra play in the preferential 
electrical response of many glasses to hydrogen ions 
with the complete or partial exclusion of the effects 
of many other ions [3] 


2. Theoretical Considerations 


Two gases separated by a membrane that is 
permeable to each will diffuse through this mem- 
brane in both directions until at equilibrium both 
sides of the partition contain mixtures of identical 
composition. Such a condition also holds true for 
a membrane separating two solutions of different 
composition, if the membrane is permeable to both 
the solute and the solvent. However, for a case 
in which the membrane is impermeable to one of 
the ions, these simple conditions no longer prevail. 
In the case of the leached surface of many silicate 


nd of this paper 


! Figures in brackets indic the literature reference t the « 


1¢ glass-solution interface 
The preferential leaching of alkali and other cations from 
s solutions, leaving a silica-rich layer to act as a nonmi- 
caused an uneven distribution of the migratable ions of the 
e dictates of the Donnan membrane equilibria and the “equation 
listribution of Ag(NH;3)%, and Br 
followed potentiometrically for Corning 015, a series of Na,O-CaQ- 
The Donnan membrane theory is compatible with 
observed phenomena such as the swelling of the 
centrations of electrolytes, 
voltage departures in the 

Donnan membrane 
glass electrode 


were investigated in connection 


ions for various experimental 


glass surface, 
voltage departures in alkaline 
superacid’’ region, and asym- 
considerations do not account for the 


to hydrogen ions with the complete or 


glasses there exists an ionogen with one ion, namely, 
the negatively charged silica-rich layer, which can- 
not migrate the glass-solution interface, 
while all the other ions of the system are free to do so 

For the simple case in which a leached glass 
surface is placed in an aqueous hydrochloric acid 
solution, the initial state of the system may lb 
represented as follows: ? 


across 


As the H* and Cl 
the interface, but the gel~ ions are not, the equi- 
librium state may be represented as 


ions are free to pass across 


To visualize the distribution of H* and Cl- ions 
at equilibrium, it is only necessary to remember 
that the oppositely charged ions must pass through 
the interface in pairs, otherwise electrostatic forces 
will be immediately set up which would prevent 
further diffusion. Therefore, the rate of migration 
of these ions across the interface depends upon the 
frequency with which they chance to arrive in 
pairs, which in turn is determined by the product 
of their concentration. Thus at equilibrium, when 


attempt will be made to clarify the confused 


?In the present derivation no 
One is merely concerned 


silicate chemistry for glasses of varied compositior 
with the facts that the silica-rich layer has a negative charge, adheres to the glas 
surface, and is permeable to the ions present in the aqueous solution. In the 
absence of more specific knowledge concerning the valence of the anions of t! 

silica-rich laver, it will merely be represented bv a single negative charge. How 


ever, any other negative valence ascribed to it will not alter the conclusior 





ransfer of H* and Cl ions from the outer solu- 
the silica-rich layer is equal to that from 
ss surface to the outer solution, the product 
concentration of these ions must have the 
value on either side of the interface 
(H*)].<[Cl-].=|H*], X[Cl-}, 
e addition of other monomonovalent electro- 
vtes to the system does not destroy the validity 
of the “equation of products.” Further, the addi- 
t polyvalent ions does not greatly complicate 
tuation. If a polyvalent ion arrives at the 
ce, it will pass across only when an equiva- 
imber of oppositely charged ions appear there 
same time. 
equilibrium, since silica gel~ as well as H* and 
ns are in the silica-rich layer, the hydrogen-ion 
tration, [H*],, in the gel must be greater than 
hloride-ion concentration, [Cl-],. The nature 
\is unequal distribution of ions may well be 
hasized and rationalized further. For the outer 
tion, let z=[H*].—[Cl-]., and for the silica-rich 
let y=[Cl-],, and lgel-],. Hence, y+z< 
As the equation of products must be fulfilled, 
y(y+-2). In this relation the product of equals 
s been equated to the product of unequals. From 
nspection it is obvious that the sum of the unequals 
greater than the sum of the equals, that is, 
2r. The validity of this is readily illustrated 
the following simple numerical example chosen 
random 
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From this it becomes evident that the concentra- 
tion of ions in the silica-rich layer at equilibrium is 
ereater than in the outer solution. This unequal 
distribution of ions gives rise to an appreciable 
osmotic pressure across the glass-solution interface. 
lf this excess concentration is represented by e, 

2r+-e. If z is eliminated by solving simul- 
neously the equations 2y+2=2r+e and 
r=y(y+z2), then z=y+ ez. This relation shows 
at a glance that z is greater than y, that is, the 
concentration of Cl- ions is greater in the outer 
solution than in the silica-rich layer. It is evident 
from this that there should be not only a difference 
in osmotic pressure, but also a difference in electrical 
potential, E=RT/F log, z/y volts, across the glass- 
solution interface. 

Without continuing further into the derivation of 
the more complicated thermodynamic aspects of the 
heterogeneous equilibria to be expected at the glass- 
solution interface, the above considerations suggest 
pertinent experiments for partially confirming or 
wholly disproving the validity of the above conclu- 
For example, the osmotic pressure predicted 
across the glass-solution interface should give rise to 

istent forces which would tend to cause swelling 
he glass surface, the extent of which would be 


S s 


partially determined by the mechanical strength of 
the silica framework and its chemical durability. 
Any such swelling or other surface alterations can be 
readily followed by means of the interferometer [4, 5] 

Another set of experiments which can be readily 
followed potentiometrically is the uneven distribution 
of Ag( NH, 3 and Br~ ions between the outer solution 
and the silica-rich surface of leached glasses 


3. Experimental Procedures 


Surface alterations of the glasses brought about by 
exposure to the various solutions were followed on 
polished specimens by means of the interferometric 
procedure previousiy used in determining the chem- 
ial durability of eptical and other glasses [4, 5] 
The exposures were made at 80° C by partial immer- 
sion of rectangular glass specimens having at least 
one surface sufficiently flat to give interference bands 
when placed under a fused-silica optical flat. This 
procedure leaves one portion of the specimen unex- 
posed for use as the reference plane. In all cases, 
swelling was plotted in the figures as negative attack 

The titration of Ag* and Br~ ions was made 
potentiometrically, using the Ag, AgBr electrode and 
the glass half-cell as the indicator and reference 
electrode, respectively. The electrode was 
adopted as the reference electrode in order to elimi- 
nate the possibility of contamination of the experi- 
mental solutions with chloride ions from the conven- 
tional salt bridge of the saturated calomel half-cell. 
Further, it served as a reliable indicator in the 
neutralization of excess NH,OH with sulfuric acid 
immediately prior to the titration of silver or 
bromide ions. The use of the glass electrode as the 
reference cell allows the potentiometric readings to 
be obtained without transference [6]. However, this 
latter feature is an unnecessary refinement, and con- 
tributes no particular scientific advantage to the 
present investigation other than operational simplic- 
ity. By using the potentiometric procedure it was 
possible to titrate the solutions in the presence of 
powdered glass samples and the AgBr precipitated 
during neutralization, thereby eliminating the use of 
color indicators and the need for filtration, washing, 
and evaporation. All potentiometric measurements 
were made at room temperature, using a Beckman 
pH meter, model G. 

The saturated solution of silver-ammonium 
mide, Ag(NH;,),Br, was prepared by dissolving sil- 
ver bromide in full-strength aqueous-ammonia solu- 
tion. Full-strength ammonia was used in every case 
because of its minimal attack on the silica-rich layer 
as indicated by previous experiments on voltage 
departures of the glass electrode [7] and the chemical 
durability of glasses [5, 8]. By taking advantage of 
the differences in solubility products of silver iodide, 
silver bromide, and silver chloride (approximate val- 
ues 10-"*, 8X 10°", and 10~"°, respectively), the silver 
bromide was precipitated and purified by the addi- 
tion of silver nitrate to a solution containing an 
excess of potassium bromide. The initial portion 
of the precipitate was discarded in order to remove 
any iodide in the potassium bromide, while gall 


glass 
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ehlorides remained in solution after 
the AgNO, and were eliminated 


the 


addition of 
The silver bromide 


uneven distribution of migratable ions at the ¢! 
solution interface for chemical Pyrex and fused si 





was washed repeatedly in hot distilled water with table 2 Thus, in the present experiments we f 
Vigorous agitation with a centrifugal stirrer to re- | surface alterations in the acid solutions ranging fr: 
move iN NO and excess Bi ions Vivorous to not detectable swelling depending 

{ll powdered glass samples were prepared by the nature of the glasses. Such results are at | 


compatible with the predictions based on Dom 
membrane equilibria. Swelling has been found 
occur for so many silicate glasses in acid buffers ¢| 
one is inclined to conclude that swelling always takos 
place in the earlier stages of attack except for so 


er ishing the experimental glasses in a steel mortar 
ind that portion was used that fell between the limits 
Standard sieves Nos, 50 and 150 


of US 


4. Results and Discussions . 
tions forming soluble silicate or other soluble sili: 


compounds [5, 8 


4.1. Surface alterations 
a 7 - , 7 = - 7 + 
a. Swelling of Glass Surfaces in Acid Solution TT 
The derivations presented in the section on theo- Na,O-PbO-Sid 
retical considerations pointed out that the concen- 
trations of ions in the siliea-rich laver at equilibrium 
are greater than in the outer solution, and that such 3 
unequal distribution must result in an osmeatic pres- 
sure across the glass-solution interface. This os- 
motic pressure must in turn tend to cause swelling 
of the glass surface in all solutions in which the . 
siliea-rich layer is formed, and the framework is s Corning ois | 
chemically durable. Table 1 and figure 1 illustrate € 2) , ' 
» » re 
typical examples of the surface alterations of several - 
; ‘ § 7 Pyrex 
glasses upon exposure for 6 hours at SO° C to a series 2 
of buffer solutions covering the range from pH2 to 2 
: = 
pH1i2. Swelling is exhibited by four of these glasses 3 
5 ; 
in the acid buffer solutions in which silicate ions are ® it cai 
not readily soluble. However, for the alkaline range s 
: = 
in which soluble silicates are formed, the rate of E ] 
solution of all the glasses exceeds the rate of swelling | Fused | 
y’ Silica 
Glasses having superior chemical durability, namely, L 
fused silica, Pyrex chemical glassware 774, and a y 
. m | >, SS eee 
Na,O-CaO-SiO, (18:10:72) glass, showed no detect- - ~ < S Y 7 
able swelling in the acid buffers during the period of - : iat es _ 
the test. Such superior glasses are known to show | 
relatively little leaching of alkali by acid solutions | 
[9], and thus the mechanism for initiating the uneven f : a ’ i 
distribution of migratable ions with the attendant ee ee ee ee ee es es ee ee ee ee 
swelling is minimized. However, it must not be ° 2 4 6 8 10 12 4 
assumed that these effects are completely absent pH 
Such an assumption would be incorrect, as will be | p : 
Ficure 1. Surface a teration-pH curves for typical si 
demonstrated by the data presented later on the glasses for exposures of 6 hours at 80° C 
TaRLe 1 Typical surface alterations exhibited ous silicate glasses upon exposure for 6 hours at 80° C to Pritton U nive 
buffe overing the range pH >to pH 11.8 
‘ Surface alteration at pH 
\ 0 kK ) Phe Cal) Ref) } ; F Pry 0.2 11.8 
Prin Fringes Fringes Fringes Fringes Fringes 
2 well 6 swell 6 swell 6 swell 1/10 swell 12D attack 
" 6 swel 6 swell 0 swell 0) swell Pitted 193 attack 
Is i swell 4 swell i swel 1/4 swell D attack 4 attack 
z2 . 0—swell 2/10—swell 2/10—swell SC e 1/4+-attack 2 attack 
sd ND ND ND Dt 1/2 attack } attack 
‘ ‘ 8 ND ND ND D 1/4—attack 1 3/4—attack 
” ND ND VNbD ND ND 1/2—attack 
eCor 0 I hen lasswa ict ! ND letectable SC, surface cut; ' D, detectable 
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b. Repression of Swelling of Corning 015 Glass in Acid Solutions 


lfeet of Concentration In the section on 
theoretical considerations the relation 2°=y(y+ 2) 
expressed the distribution of ions on both sides of 
the interface The potential difference across the 


interface resulting from this uneven distribution 
was given as E=RT/F log, x/y volts This relation 
is true even if other ionogens are added to thesystem, 
including neutral salts or other electrolytes having 
polyvalent ions [2]. The addition of anv electrolyte 
must necessarily produce a change in the potential 
difference by disturbing the equilibrium, yet all ions 
present when equilibrium is again established are 
producing the same potential difference regardless 
of valence. This latter condition may be difficult 
to visualize, but appears directly in the derivation 
involving polyvalent ions [2] (see p. 435 on the 
effect of valence 
If the equation r=yly >) is solved for y and then 
substituted in E=RT/F log, x/y, the relation for the 
potential between the silica-rich layer and the outer 
solution becomes 
») 
E=RT/F log, = 


y 4a 


d 








From this it is obvious that upon increasing the 
value of ¢ while 2 (the negatively charged silica lave1 
remains almost constant, / must decrease (along with 
the osmotic pressure and resulting swelling), ap- 
proaching zero as a limit, that is 


») 
J 


limit K=RT/F log ‘ 0 


To test the validity of this reasoning as applied 
to glasses, Corning 015 electrode glass which ex- 
hibited a moderate rate of swelling, table | and 
figure 1, was exposed to various concentrations of 
organic and mineral ne ids la order to emphasize 
the effects, the exposure time was tacreased to 48 
hours at SO° C The results ure plotted in figure 2 
In every case, the rate of swelling was decreased by 
the higher concentrations and approached zero for 
Aces rding 
to the theory for swelling as applied to gelatin, the 
swelling should go through a maximum in passing 
from vers dilute to higher concentrations of elee- 
trolvtes 2,10} This same reasoning also holds true 
for silicate glasses. However, earlier experiment 
with the swelling of Corning 015 had not revealed 
such a maximum for a concentration series 10 \, 


107? N, 107' N, and N [5] 


the concentrations above 50° percent 


2 Effect of vale nee upon swelling Although the 
data plotted in figure 2, A, were obtained on solutions 
prepared only on an approximate percentage basis, 
they indicate that the acids differ in their capacity 
A more carefully 
controlled set of experiments was carried out with 
solutions prepared on a molar basis The results 
are plotted in figure 2, B. The values for swelling 
brought about by the monobasic mineral acid, 
HCl, and the monobasic organic acid, CH,COOH 
fall essentially on the same curve, whereas the 
swelling in the case of the dibasic acid, H.SO,, is 
appreciably less The effect of the tribasic acid 
H,PO,, follows in a general way the monobasic acids 
All of these results can be rationalized, provided it 
is assumed that the silica-rich laver is not dissolved 


to repress the swelling of glass 


by any of the acids 

If the heterogeneous equilibrium at a glass- 
solution interface is a characteristic Donnan mem- 
brane equilibrium established by a nondiffusible 
negative charged gelatinous silica, the repression of 
swelling of the glass surface brought about by 
monobasic acids, such as HCl and CH,COOH 
should be essentially the same. On the other hand, 
strong dibasic acids should repress the swelling more 
strongly The present procedure is not sufficiently 
exact to give a quantitative verification of these 
predictions, but the results are generally in accord 
with them The tribasic acid, H PO,, dissociates 
as a monobasic acid below pH 4.7 [11] and has been 
shown in the case of gelatin to influence the mem- 
Lrane potential as well as the osmotic pressure and 
swelling the same as the monobasic acids 2] 


From these observations it is evident that the commercial practice toring 


concentrated acids in glass cont 
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4.2. Uneven distribution of Ag(NH,)> and Br 
Ions Between the Silica-Rich Layer of the Glass 
Surface and the Outer Solution 


In the previous sections the uneven distribution of 
ions at the glass-solution interface has been assumed 
demonstrated As neither the 
glass surface nor the 


but not actually 
concentration of ions within the 
differential existing across the 
readily obtainable by 
to design experiments that 


voltage vlass-solution 


direct: measure- 


interlace 1s 
ment, it Was 
demonstrate the excess of the 


hecessary 
would migratable 
positive ions appearing in the 
or the excess negative appearing in 
phase Ag(NH sr dissolved in full-strength aque- 
NH satisfactory indicator 
for such experiments. Full-strength NH 
because of its low attack upon the silica-rich laver 
although the pH is high [5,7] 

Glass samples (approximately 10 ¢ 
a U.S. Standard 50-mesh sieve and were retaimed on 


swollen glass surface 
ions the outel 


Ons 


ous provided very 


was used 


whic h passed 


a 150-mesh screen were leached at pH 2 tor 6 hours 
at 80°C, with occasional stirrin Thev were then 
decanted and rinsed twice in full-strength aqueous- 
ammonia solution, followed by 5 minutes in 50 ml 
of saturated Ag(NH,).Br solution with intermittent 
agitation The glass samples were then separated 
from the liquid phase by decantation 
and liquid phases were acidified with concentrated 
H.SO, and titrated potentiometrically for 
Ag* remaining in the gelatinous silica-rich surface 
and the excess Bi the outer 
solution For solutions 
were made strongly acid to insure that little change 
in pH takes place during the titration, and also to 
ready coagulation of the AgBr pre- 


The glass 
excess 


rons 
determinations the 


appearing in 


these 


bring about 


cipitated. 

a. Na,O-CaO-SiO, Glasses 
titrations were made of the 
Ag(NH,)? ions appearing in the glass surface in 
excess of the Br ions for a series of Na.O-Cal -Si0, 
glasses containing 70 percent of SiO, and various 


Potentiometric 











is rs at si° ( 


percentages of Na,O and CaO. The data are giv 
in table and summarized in figure 3, where 
swelling of the glasses at pH2, and the ex: 
\g NH ions appearing at the class surface 
plotted against the composition (percentage of CaO 
of the The results are 
of the low titration values obtained on the classes 


classes interesting because 
superior durability as compared with the glasses 
inferior durability. Another interesting feature 
these data is the fact that for the two glasses 
maximum swelling (29.5 and 29.75 percent of Na.O 
respectively), the excess silver appearing in the g 
surface decreased, falling essentially to zero for tl 
end member of the series, 30 percent Na,O This 
apparent anomaly will be discussed and clarified | 
the data obtained on the Na,O-SiO 


aa b 


classes sectiol 


b. Na,O-SiO, Glasses 


A similar set of data were obtained for a series 
soda-silica glasses, table 2 and figure 4 The excess 
Ag NH ions appearing in the negatively charg: 
silica-rich phase increased for the glass compositions 
showing increased swelling down to the composition 
of 74 percent of SiO,. Below this percentage 
silica, the Ag NH ions retained at the interfac 
fell sharply toward zero, as in the Na,O-Ca0-SiO 
serics, because the swollen silica-rich surface of the 


specimens was decreased by solution and by sloughing 


off of the mechanically weak gel. These data illus 
trate well the difficulty encountered in obtaining con 
sistent durability data by the interferometer procé 
dure in studying the behavior of very poor glasses 
acid solutions. Specimens of these glasses given 
minute exposures suffered attack as shown in figui 
4, whereas exposure of the same glasses for 5 minutes 
or less resulted in liberal swelling [12, 13]. 


c. Effect of pH 


The data for the excess Ag(NH;)? appearing 
the glass-solution interface after samples of Corni: 
015 were leached in buffer solutions of different p 
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Chemical Durability pH 
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Excess Ag (NH 
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Chemical Durability 
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ristobolite 


Chemical Durability pH z 
mirumurm exposure 


listed in table 3. The resulting curve * is com- 
pared in figure 5 with the silicic-acid titration curve 
1], the chemiecal-durability 5] and voltage-depar- 
re [3, 15] curves of this glass over the same pH 
\ set of chemical-durability data that is rep- 
resentative of the “‘superacid”’ region is included. It 
s obtained with sulfuric-acid solutions and reported 
log 1/[H.SO,]. From inspection of these curves it 
vident that as the negative ly charged siliea-rich 
iver Is carried away as soluble silicates by alkaline 
solutions, the magnitude of the uneven distribution 


range 





40 , aa 20 
36} 
- 
Oo 32} 
- 
o 26 3 
S 27 2 
9 Excess Ag(NH_J* m = 
‘ 2 0 u 
= <V . 7+ < 
+-™ 6 9 408 o 
> = 
x = 
z 2 . 06 < 
o ‘ 
a 6 ,04 
4 4 we Tregted After " Jo2 
* Leoching of pH6 
0} , 40.0 
?) Y 
30 , { O loz 
s oo Log/(n,SO,) Cremco! Durability a 
$ 0 } 
& 
° a Voltage Departure 
ry 0 +0 
5 2 Ne,0 S10. 
2 a P S10, 
8 5 
& -20} map 210, | 40 | 
‘ 
na) 
30 ’ 
Na,0 25:0. 4) 
-40} a ® is0 6ClO 
SiO, ° 
Free #Ci Titration, Silicle Acid E 
120 
| Jieo 
2 ie) 2 4 6 8 2] 2 4 
pH 
Fiaure 5. Effect of the pH of the buff ) ms on the chen 
dural voltage dep e, and the excess Ag(NH 
pared with 


ons appea ng n the giaas 0 or ’ } omi 


of migratable ions at the interface decreases This 
condition readily rationalizes the controversial, so- 
ealled “‘nonreversibility of the glass electrode [8, 
16, 17, 18}. The removal or partial removal of the 
silica-rich laver by alkaline solutions must necessarily 
cause a voltage shift as the nonmigratable member 
of the interfacial system | changed The siliea-rich 
surface is not whollv regenerated until the electrode 
is immersed again in an acid solution 


ras.e 3 Effect of the pH of the leaching solution upon the 
surface alteration and upon the uneven distribution of migrat 
thle ons appear:ng at the surface of Corning OLS glass 
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In the low pH range, the correspondence of the 
repression of swelling and the uneven distribution of 





migratable ions in H,SO, solutions (table 4 and fi 
with the voltage departures reported for the “s 


acid” region, is entirely logical. However, ther 
several other considerations that must not be ¢ 
looked First, it must be borne in mind that 
decrease in the water activity of the solutions at 
increased concentrations also parallels the perf 
ance ol the class electrode and strongly supports 
theory that the glass-electrode functions as a “‘w 
electrode 3, 19, 20]. This theory is also stror 
supported by the general correspondence of the 
groscopicity of experimental glasses with their 
trode function (pH response) [13, 21] 

However, as pointed out in the present invest 
tion, there should be a lowering of the interfa 
voltage with increased acid concentration in ac 
with the equation 


limit k—RT/F log _ 0 
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This voltage shift should theoretically be independe: 
of any water activity interpretation, but obvious 
for aqueous solutions it would be accompanied by 


decrease in the water activity of the solutions, an: 


the two theories would have this same featur 
common Yet the source of the voltage depart 


in the first theory would be internal, whereas in th 
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im the “su pe rac id”’ region and the accompanying unever 
tribution of migratable ions 














of 


second, the postulated cause of the voltage departure 


would be at the interface 








( y OT gli spon erpo 7) tS hours at SO C’ to 
nt ms of aq H.SO,, and th never . 
qd i , ng } gd 

} H 

HI 
\ further consideration, the repression of swelling 
vith tl iccompanving decrease in the formation of 
migratable ions and the attendant displace ment 


the heterogeneous equilibria at 
essentially 
ost ated DY the 


the glass-solution 
a byproduct of the conditions 
equation 
In selecting a theorv that ade quately accounts for 
the 
region there is a liberal list from which 
namely, (1) Dehydration of the 
ass 8, IS] 2) alteration of the physical and 


i voltage departures of the class electrode in 


| 
superacid 


o choose ‘lectrode 


chemical nature of the surface of the electrode [5. 8] 
}) operation of the glass as a “water electrode > 
q ”!) and (4) alteration of the interfacial volt- 
g cord with the above equation [1, 2 If any 

of them ts correct, any disagreement as to which 
s can merely center about agreement on which 


mm ver of the list is the cause and which ones are 
Probably. as sso prevalent throughout 
anv ultimate explanation will be faced with 


al nomaly or contradiction 


the effeets 


nature 


4.3. Surface Equilibria and Glass Electrode Voltage 
Alt ho igh 


such a 


Donnan membrane phenomena occupy 
role in the glass-solution interface 
equilibria, they apparently can play little or no part 
in the preferential response of the 
the hydrogen-ion activity of aqueous solutions and 
the partial or complete exclusion of the effects of 
the other This conclusion 
can be reached by an inspection of the basic equa- 
tions of the equilibrium state 

The equation for the voltage appearing across the 
zlass-solution interface due to the uneven distribution 
of a migratable monovalent ion is E=RT/F log, a y 
volts. If an ionogen yielding an ion of any other 
aleney, MM’ is added to this system, the voltage 

ntributed by this ion at equilibrium will be 





dramatu 


glass ele ctrode to 


ions in the system [235] 


E=RT/nk 


log 








From the equation of products \V/ q Cri 
M"*}en X (CI and) (H*|: <(Cl-]: [Hf 
C |} from which it follows that 


Therefore hk R17 I log rey 


It is evident from this derivation that at equilibri 
um all ions in the svstem 
identical 


Thus 


tvpe ot heterogeneous equilibria is concerned 


regardl ss of the il valence c 


are contributing 
interlace 


voltages across the class 
solution as far as this particular 
the 
well as an 


glass electrode should respond equally 


ton in 


the aqueous solution [24, 22], 
provided all other positive tons were held constant 
except the consideration (As this latter 
conclusion is so COnSpl¢ uously contrary to the rene 


ally 


trode as used, it 


indicator of any 


one under 


observed voltage perlormance of the lass elec 


is obviously lo look else- 


hecessary 


where for an explanation of the selective voltage 


response oft thre electrode to hvdrogen mon 


glass 
This does hot imply that the voltage ar ross 
to the 


ions is negligible or 


activils 
uneven distri 
that it 


the glass-solution interface duc 
bution of the 


plays no part in the observed voltage of the 


migratabl 
rlass 
electrode 

interfacial voltage seems to be of 
For example, consider the 
that 
surtace 1s partially 0 


In reality this 
considerable magnitude 
shift 


when the 


lara voltage so-called “soda error 


silica-rich 


appears 
completely removed by alkaline 
20) that form 
obtained on the voltage 


solutions > ) 19 
From similar evi 
shift in hvdrofluort 


acid solutions, it appears that the interfacial voltage 


soluble silicates 


dence 


caused by the presence of the gelatinous silica-rich 
laver must often be than 250 my 
These data with a 
glass electrode of ( ornming O15 compared augcainst the 
hvdrogen electrode [3 p. 124; 25 The 
voltage de parture cannot be attributed to hvdrogen 


figure 7 
well-conditioned 


create 


were obtained 


obser \ ed 


1oOns, because both electrodes when unattacked re 
spond equally well to hvdrogen-ion activity Fur- 
thermore, the voltage shift cannot be attributed to 


an equilibrative either electrode to 
fluoride neutral no such 
response to fluoride ions is observed [3, p. 131 The 
complete or partial removal of the silica-rich layer 
with its interfacial distribution of migratabl 
apparently the main contributing variable 


response ol 


ions. because in solutions 


Ons Is 


The surface equilibria also play an important role 
in the origin of the asymmetry pot ntial of the 
The miner 


glass 
and outer surface of a glass- 
bulb of the Cremer-Habet tvpe 25 26) 
are known Lo be different, as shown by their 
capacity With different hetero- 
geneous equilibria and accompany ing voltages of 
different magnitudes on the two surfaces an asym- 
metry potential for the glass electrode is inevitable 
Further, upon taking appropriate steps to make the 
surfaces identical, the asymmetry potential 


electrode 
electrode 
very 


to sorb dves [17] 


two 
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Preparation of )-Mannitol-C ‘ and Its Conversion to 
)-Fructose-l-(and 6)-C'' byAcetobacter Suboxydans 


Horace S. Isbell and J. V. Karabinos 


Methods are described production of p-mannitol-1-C™ and p-fructose-1l-(and 
6)-C™ Thess ateria 1 wowerful new tools for research in biology and chemistry 
\la tol-1-¢ Vas prepa S0-percent vield from p-mannono-y-lactonm p-Fructose- 
l-(and 6)-C™ was produced 4-percent vield by growth of Acetobacter suboxydans on a 
\iar }-( ! 1 
l. Introduction As both p-mannitol-1-C™ and p-fructos 
and 6)-C™ would provide powerful tools for tl 
The Bureau recently undertook a program for the elucidation of many problems in biology and chen 
istry, the preparation of these compounds w 


development of methods for the production of a 


variety of position-labeled sugars || The prepara- 
tion of p-glucose-1-C™ and p-mannose-1-C™ in high 
Vie ld was cdlese ribed In a previous re port 2| This 


D mannitol 1-C"* 
and 6)-C™ by 
suboxvdans Priot to 
labeled fructose 


paper deals with the produc tion o 
eonversion to pb fructose ] 
Acetobacter 


randomly 


and its 
oxidation with 


the present investigation 


had been prepared by photosvnthes S but posi- 
tion Leaboe led fructose had not been reported It was 
known that p-mannitol could be obtained by redue- 


and that it could be converted to 


\c ‘ tobac ter suboxvdans 


tion ot D-mannose 


p-fructose by oxidation with 


1] p-Mannitol differs from p-mannose in that the 
ends of the molecule are alike For this reason a 
l labeled mannitol | can be considered is either 


l- or 6-labeled, and anv unsvmmetrical derivative 


prepared from this substance would be 1,6-labeled 
Thus 


provided the carbon chain remains intact 
obtain 


by oxidation of p-mannitol-1-C'™ one would 
and 6)-C"™ (II 


D frie tose-|- 


H H H OH OH i 
Hor ( ( ( ( ( OH | 


H OH OH H H H 


6 ) i } 2 
Mar 1-( 
H H H OH oO H 
Ho ( ( ( ( C*—OH and 
H OH OH H H 
H H H OH O H 
HoO—c* ( ( ( ( ( OH I] 
H OH OH H H 
p-fructose-l-(and 6) -( 
r At Ey 
, {the pat 
en ed 
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undertaken 


2. Discussion of the Processes 


includes 


The production of p-mannitol-1-C"™ 
by tl 


p-mannono-y-lactone-1-C"™ 
described [2], (2) conversion of tl 
hy sodium amalgan 
reduction catalytic reduction of the suga 
with hvdrogen. Each step in the 

studied in detail, and procedures were developed 
p-mannitol-1-C™ from p-mannono-y-lacton 
all radiochemical vield of about 80 


preparation of 
method already 
lactone to p-mannose-1-C"™ 
and (3 
process Was 


obtain 
1-C* in an over 
percent 

The production of p-fructose-1l-(and 6)-C"™ involves 
oxidation of p-mannitol-1-C™ with Acetobacter 
oxvdans It has been established 5] that Aceto 
bacter suboxvdans preferentially oxidizes the 


sub 


OH OH 
R—¢ ( CH,OH 
H H 
group to the 
OH () 
R—¢ ( CH,OH 
H 


group. In agreement with this property, p-fructos« 
l-(and 6)-C™ is obtained by oxidation of p-mannitol 
with Acetobacter suboxydans.‘ As pointed out 
section 1, the ends of the molecule are alike. Fo 
this reason, oxidation could proceed equally wi 
from either end of the molecule, and, in the ever 
that oxidation took place at both ends, a diket: 
derivative would be formed. Thus.  prolonge« 
action of Acetobacter suboxydans on mannito 
should vield the diketo derivative 


H OH O 
CH,.OH—¢ ( Cc C—CH,OH 
O OH H 


p-threo-5-0x0-2-ketohexose 


4‘ See footnote 




















fact 
a series of man- 


To ascertain whether the oxidation did, in 
proceed bevond the fructos 
tel Acetobacter suboxydans cultures were set up 
ind. as given in table 1, the 
was followed by determinations of fructose and total 
copper-reducing substances The fructose 
ration was determined from the polarization at two 

mperatures 6] 

rmined by the modified Scales method 7| 


of figut 


stage 
course of the reaction 
concen- 


Wiis de- 
Curve 


Total reducing substance 


y | gives the vield of p-fructose determined 
from the difference in optical rotation at two temper- 
at Curve 2 gives the copper-reducing power of 
solution expressed as the percentage of the re- 

power of a solution containing the amount of 
ictose molecularly the mannitol 
originally present. The copper-reducing value repre- 
sents all reducing substances present and is expressed 
The vield of p-fruc 
tose ren hed a maximum in about and then 
decreased fol about 6 davs when active growth of the 
sm had slowed down In every sample, after 

he first day the copper-reducing value exceeded that 
of the p-fructose determined from the polarization ut 
wo temperatures Hence, the 
nt having copper-reducing powet other than 


fy equivalent to 


arbitraruyv on the fru LOSe S« ale 


2 days, 


solution contained a 


onstitu 


»-fructose Presumably, this substance was the 
imticipated diketo compound or a substance de- 
ved therefrom The results clearly show that it is 
ssarv to control th reaction time to «6obtain 
im vields of p-fructos The experiment of 


th a surface-volume ratio 
egreement with the work of Fulmer, 
Dunning, and Underkofler, a somewhat higher vield 
f f tose was obtained with a higher surface-volume 


| ! was conducted w 


ol Ss In 


In the preparation of p-fructose-1-(and 6-)C" 
tracer level, a surface-volume ratio of 1.59 was 
ised, and a radiochemical vield of 54 percent was 


' | 
ed] 


3. Experimental Procedures 
3.1. D-Mannono-y-lactone-1-C'' 


This substance was prepared in 67-percent vield 
by t method 


sodium-bic arbonate Cul bon-cdioxide 


»)] 
| 


pre \ ously deset ibed 


3.2. Conversion of D-Mannono-y-lactone-l-C to 


D-Mannitol-1-C" 


One millimole of p-mannono-y-lactone-1-C'* having 
an activity of 10 we was dissolved in 20 ml of water; 
> g of benzoic acid was added, and the mixture was 
stirred in. an ice bath until cold. Sodium amalgam 
pellets (4.6 ¢ contaming 5 percent of sodium) were 
dded and stirring was continued. After 2 hours the 
rcury was removed, and another addition of 1.2 
benzoic acid and 


g 
1.6 ¢ of amalgam was made 
rring was continued at 0° to 5° C for 2 hours, after 
vhich the mercury was removed, and 1.26 g of oxalic 
acid was added. Then the mixture was extracted 
vith chloroform to remove the benzoic acid The 
sodium oxalate in the aqueous solution was precipi- 


wy 








et 
1 
| 
ACT TET 
hy RE | ij ? ni | } 


tated by the addition of approximately 2 volumes of 
methanol and 2 of ethanol The insoluble salts were 
separated by filtration and washed with methanol 
The resultant aleoholi liquor was concentrated under 
reduced pressure to about 5 ml and was diluted again 
with 2 volumes of methanol and 2 of ethanol The 
resultant precipitate was separated and washed with 
methanol 
apy roximately 


The alcoholic liquor was evaporated to 
2 ml and was then transferred with 3 


ml of water to the tube of a high-pressure hydro 
genator containing 0.1 ¢ of calcium carbonate and 
0.5 ¢ of Raney nickel catalyst The mixture was 
hvdrogenated for 2 hours at 125° C and 1,000 Ib/in 


Then the mixture was cooled to room temperature 
filtered, and passed through a column (1.4 by 25 em 
contaiming equal portions of Amberlite IR 100-H and 
Duolite A4. The effluent and wash liquor were com- 
bined and lyophilized. The residue was dissolved in 
the minimum quantity of hot methanol After the 
addition of an equal volume of isopropanol, and 
seeding, the solution was set aside for crystallization 

In the course of 3 davs 117 mg of 
"4(6.4 we It melted at 
in substantial agreement with the accepted 

166° to 168° © “Carrier”? p-mannitol 


of p-mannitol 
p-mannitol-1]-( 
165 he 

value of 


was obtained 


100 mg) was dissolved in the mother liquor and 
after crystallization, 98 mg of p-mannitol-1-C'* was 
obtained with a radioactive content of 1.6 ue. The 


total radioactive recovery was S80 pereent 


3.3. Preliminary Study of the Oxidation of D-Mannitol 
to D-Fructose with Acetobacter Suboxydans 


Aliquots (25 ml) of an aqueous solution containing 
9.1 g of p-mannitol, 2.5 g of veast extract, and 1.5 g 
of potassium dihydrogen phosphate in 500 ml were 
placed in 125-ml Erlenmeyer flasks and sterilized for 
15 minutes at a steam pressure of 15 1b. Each flask 
was then inoculated with five drops of a previously 
prepared inoculum of Acetobacter suboxydans and 
allowed to stand at 30° C in an incubator. The 
inoculum was prepared by growth, at 30° C for 45 
hours, of a fresh culture of Acetobacte suboxvdans 
on a solution containing 2 g of mannitol, 0.2 ¢ of 
yeast extract, and 40 ml of water 
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following the inoculation 
from the incubator 
25-ml volumetric 
One-half 
then 
centrifuge 


At regular intervals 
one of the removed 
The liquid was transferred to a 
flask diluted to volume 
neutral lead acetate trihvdrat 
Th 


and the 


flasks was 
gram of 
added 
tube. 
centrifugation 


and 
was 
transferred to a 
was removed by 


mixture Was 
precipitate 
To the clear solution was added 0.25 ¢ of potassium 
and the resulting precipitate 
The clear 
for measurements of 


oxalate monohydrate 


was again removed by centmfugation 
solution was then used directly 


opti al rotation and coppel -reduc Ing powe! 


PAR. 


The optical rotations given in table 1 were made 
at two temperatures with a Bates saccharimeter and 
a 2-im tube. The rotations at the higher tempera- 
tures were corrected [6] for expansion by use of the 
coefficient 0.00044 for each degree rise in temperature 
The difference between the corrected rotation and 
that at 20°, divided by the temperature difference 
and by 0.03441, equals the number of grams of pb- 
fructose per 100 ml of solution. Complete con- 
version to p-fructose of the mannitol originally in the 
solution would have given 1.8 g of p-fructose per 100 
ml 

Copper-reducing values, also given in table 1, were 
determined on the solutions by the modified Scales 
method [7] standardized against p-fructose. The 
values so obtained represent all substances that re- 
duce the alkaline copper reagent 


3.4. Conversion of D-Mannitol-l-C'' to D-Fructose-l- 
(and 6)-C"* 

One millimole (182 mg) of p-mannitol-1-C'*, with 
a radioactive content of 8.8 yue/millimole, was dis- 
solved in 10 ml of aqueous solution containing 50 mg 
of yeast extract and 30 mg of potassium dihydrogen 
phosphate in a 50-ml Erlenmeyer flask (surface/vol- 
ume= 1.59 cm?/ml The mixture was allowed to in- 
cubate at 30° C with three drops of a previously pre- 
pared inoculum of Acetobactor suboxyvdans to the 
point of maximum rotation and reducing power (45 
hours). To this solution was then added 1 ml of a 
20-percent solution of zine sulfate and an equivalent 
quantity (previously determined with bromothy- 
mol blue indicator) of a saturated solution of barium 


hvdroxide. After removal of the precipitate by « 
trifugation, the solution and washings were pas 
slowly through a (1.2 by 22 em) column contain 
a mixture of Amberlite IR 100—-H and Duolit: 
The effluent and aqueous wash liquor were con 
trated by lvophilization to an amorphous resi: 
weighing 145 mg. This residue 
several drops of methanol, isopropanol was ad 
to the point of incipient turbidity, and the mixt 
was seeded with p-fructose. Crystallization was 
lowed to proceed for several days, after which 56 
of crystalline product was obtained. The spe¢ 
rotation ({al]? 82°) as well as the radioactiy 
0.0428 ue/mg) of the material was about 10 per 
too low Hence, the \ ield of C ‘labeled fructose wv 
actually 49 mg (27%) instead of 56 mg. An ad 
tional quantity of labeled fructose was separated 
the use of 200 mg of carrier p-fructose. An aque 
solution of the carrier was added to the mother liquor 
and the solution was evaporated to substantial di) 
The residue was taken up in 1 ml of methane: 
and crystallization was induced by the addition of 2 
ml of isopropanol and seeding with crystalline fru 
The crystals that separated weighed 210 mg 
and had an activity of 0.0114 ue/mg The crop of 
fructose undiluted with carrier contained 2.40 we, and 
the carrier-diluted sugar contained 2.39 ue: thus the 
radiochemical vield was 4.79 ue, or 54.4 percent 


was dissolved 


hess 


tose 


The authors express their appreciation to Nancy 
B. Holt and to Abraham Schwebel of the Bureau, 
who conducted the analytical work, and to Laura C 


Stewart of the National Institutes of Health, who 
supplied a culture of Acetobactor suboxvdans and 


made suggestions concerning its use 
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Analysis of Goniophotometric Reflection Curves 


Isadore Nimeroff 


appearance apart trom 


ar distribution of the refle< 
shows that ess¢ 
hod of 

pretation 
Phe 
of the 

dex is 


Statistics 
whereby sucl 

dex is show! 

images 


KUTLOSIS In 


reflected 


1. Introduction 


rsuit of a method for specifying glossiness, 
ster, of surfaces many theorists, researchers, and 
ieers, physical and psychophysical, have devoted 
oht and time to theorizing and experimenting 
hniques for measurement have been devised, 
ning irom observation with a single set to obser- 
many of illuminating and 
viewing conditions. The former of these techniques 
oversimplified solution to the problem; the 

r is too costly and time consuming to perform 
every surface under consideration. Such a 

lete set of data was published by Thaler [1] ! for 

A compromise 
en these extremes is required. Such a com- 
se involves obtaining goniophotometric curves 
the angle of specular reflection in the plane of 
surement for one unidirectional incidence angle. 


t sets 


ions involving 


s of magnesium oxide surfaces 


2. Goniophotometers 


The earliest known goniophotometer is the one 
devised by Bouguer [2] in attempting to corroborate 
Lambert’s law [3] by experimental evidence. Since 
then many goniophotometers have been developed 
and The problem of how to analyze and 
interpret the curves obtained on goniophotometers 
remains This paper proposes using an estab- 
lished statistical method analyze monoplane 
goniophotometric reflection curves, yielding numerics 
that correlate with specular gloss and estimates of 
distinctness-of-image gloss 

Many workers in the field of glossimetry have 
realized that the only way to specify completely the 
optical properties of reflecting surfaces is by analyz- 
Ing goniophotometric curves of the reflected light 
In the interest of economy one goniophotometric 


ces ribed 


still 


to 


reflec tion curve for one angle of incidence and one 
plane of view may provide an adequate account of 
what happens to the flux reflected by a surface 
his flux distribution can be measured visually or 
hotoelectrically The most detailed description of 
i visual goniophotometer is given by MeNicholas [4]. 
His instrument had provision for diffuse and unidi- 


ions and tert 


texture, 
ted light, but a sun 
ntial parts of 
applied to th 
adata ar 


appli ad to om pre 


is defined 


negiect 
tri 


hese fundament beer 


norome 


probien 


piel 


view He 
distribution 


multiplane 
obtain 


illumination and 


goniophotome ter 


rectional 
his 
curves for diffusing media 

Wetlaufer and Scott 15] deseribed a photoelectric 
goniophotometer with rectangular source and receiver 
by 5.5 but the relative 
SIZes these apertures were riven They 
obtained goniophotometric data near the angle of 
specular reflection with 45° incidence for 13 materials 
exhibiting reflection ranging from nearly completely 
completely diffuse, and drew 
conclusions about the relative merits of large and 
small incidence angles in the measurement of specular 
gloss. Only one incidence angle was used, and no 
analysis of the goniophotometric data was made 

In a paper on 60° specular-gloss measurement [6] 
Hammond and Nimeroff described the optics of a 
phétoelectronic monoplane goniophotometer having 
circular source, and 0.6° circular 
With this instrument fractional 
were obtained for a of 
specimens differing in appearance, yet 
having approximately the same 60° specular-gloss 
readings according to the ASTM Method [7] The 
curves indicated that these specimens had radically 
different fractional-reflectance distributions If the 
information contained in curves could be 
adequately expressed as a single number, a better 
index of gloss characteristics would be obtained than 
is possible from the simple measurement of specular 


used to 


apertures summing to 50’ 


of not 


specular to nearly 


high resolution, 0.3 
receiver apertures 
reflectance curves 


series 


obviously 


these 


loss 


3. The Concept of Moments 


A method for the description of these fractional- 
reflectance curves may be found in the field of sta- 
tistics, where the technique of characterizing various 
types of frequency distributions is of fundamental 
importance In statistics the measures of dispersion, 
skewness, and peakedness have been found to facili- 
tate comparisons between distributions. These 
measures are computed by the use of moments. A 
brief résumé of the concept of moments follows 

With observed frequencies, Ii, Jo, l, 
and with 

f, equal 


responding to variables, 2, 2 
the sum of the frequencies, f, 


W 





kth moment ts 


here WN pS In other words the 


he arithmetic mean of the &th powers of the variate 
When the taken about the 


centroid of mean T or the then 


axis of moments ts 
, 


variable 


kth monent about the centroid, or mean 


3 the 


4. Curve Analysis 


Harrison |S] suggests that polar-intensity distribu 
tion curves should be analyzed to obtain the desired 
information about surface He suggests 
also that these polar curves should not be considered 
and straight lines 

spec ular components lines 
phvsically from the diffus 
Barkas [9] suggests a method for theoreti- 
fairly 


appearance 
as composed ol elliy ses, circles 
because the cannot be 
separated components 
circles 
cally 
matte surfaces by substituting an equivalent surface 
for a physical surface. The Barkas treatment has 
not vet been applied to semucloss\ and ole Ssv surtaées 
10] found that by summing the diffuse 


separating these two components for 


Pokrowski 
reflection and the specular-reflection components one 
could obtain for illumination the 
distribution of reflected flux for all directions The 
made on the 
Lambert's 


a given angle of 
computation of these components is 
assumption that diffuse reflection obevs 
law and that the specular reflection obeys Fresnel’s 
law [11] for an equal distribution of mirror facets 
oriented at all possible angles. Chinmayanandam 
[12] derived an exponential expression to account 
theoretically for the reflection from optically rough 
assuming that the angles between the 
reflecting elements and the gross plane are randomly 
distributed according to the probability law 

Typical gomophotometric c observed for 
materials ranging from matte to high gloss are shown 
in figure | Clues as to how to analvze these data 
are difficult to obtain from a polar graph. However, 
if rectangular coordinates are used for the data in 
figure 1 and fractional reflectance, F, is plotted on a 
logarithmic ordinate against angular 
departures from the specular angle, d, on a linear 
their resemblance to nor- 
about the arithmetic 


surfaces by 


irves 


scale as the 


see fig 


distributions 


scale as abscissa 
mal frequency 
mean becomes apparent 


In the analysis of frequency distributions, mathe- 
matical statistics (see, for example, Reitz [13] and 
Croxton and Cowden [14]) has shown that three 
indexes are important. These three indexes are 











ophotomet 


coordaitnawes 











d, ““2REES 


Liophotomet 


oordinates 


dispersion index, o; skewness index, aj; and kurtosis 
index, a, For the purposes of this paper, they are 


defined as follows 


The first index, ¢, gives an indication of the spread 
of the distribution, the second, ay, gives an imdicatio! 
of the dissymmetry of the distribution about th 
mean, and the third, ay, gives an indication of th 
peakedness of the distribution. This last index ts 
of primary importance in goniophotometric curv: 
analysis for, as will be shown later, its correlatio! 
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| vith estimates of distinctness-of-reflected-image computation is further simplified, because we can 
| oss is almost pertect readily take 
; \ me thod was devised by Pelton 15} inh whi h he fal { e 
j ns a series Of measurements of photometric 
ehiness, By, by tilting a specimen about an axis | where d, is the departure from the specular ang 
i the photometri plan Thus he obtained a typ assumed mean), and is the goniophotometrically 
oniophotometric distribution curve with simul- | obtained fractional reflectanes distribution — fre 
aneously varving incidence and viewing angles. | quen 
} 1} irve suggested to him an analogy to the fre- PaABLeE 2. ¢ 
que distribution curve according to the law of 
1) bility ae then computed the dispersion 
isually called standard deviation using 
yt) 
Ss, BO 
\ SB 6 
vl 0 is the angle of tilt of the specimen This 
is analogous to taking the square root of 
for which 4 2 As the standard deviation 
nverse relation to the sharpness of his curves 
on took 1l/e¢ as a measure of what he called 3 
4.1. Computation Short Cuts zs st 
lo shorten the computation we define deviates . 
rom an assumed mean Then moments x 
van | ; 
\ , 
are related to x, about the mean bv the 
expression obtained from the xpansion of 
’ nomial in eq 2 


o 4° 
m=>(, ) 8 
' 4.2. Characteristics of the Indexes 
| 7.9 . . 
( ) . , If a distribution is symmetrical, the skewness 
; Ale! s\t > , 
index, ay, has a value of zero There is no uppel 
limit to the magnitude of a Mathematically 


an integer having values, 0<t</ stated 


hi st five moments about the mean (2, m, m., 
tT and wT are im tabl | Dispersion skewness. 0) a 


and kurtosis indexes may then be defined in terms 
of these moments Values for a, as great as 2 dicate marked skew 
hess, howevet! 

The kurtosis index, a,, has a value of 3 for the 
normal curve, y=ae~? For curves more peaked 
than the normal curve, a, is greater than 3, whereas 
for curves flatter than the normal curve, a, is less 
than 3. There is no upper limit for a, For a flat 
curve, that is for f, constant, the value of ay, desig 
nated a,, Is determined from the number of obser 
vations on which computations are mace Krom 


. u the definitions of ay, F and mw. we see that 
lo facilitate the computation of the three indexes 


data were taken at the angle of specular reflection 


dat 1° intervals for 5° on either side of the specular » Dai hidt/N 24} 
ingle. To show the ease of computation an example a - ! 10) 
7 shown in table 2 for a plaque of white ceramic tile . ( S*. f.d2/N ) ( 3‘ ) 
luminated at 60°”, sper ular gloss of about 90). The —— 2 
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where 7 V7 number ol angle settings Phe 
general case for 1° intervals is shown as follows 
We ha 
S\ j >'S \ ) 
> 1! 
St Hi? S vy 
vhere A, represents positive numbers rang from 
to .V/. since observations are made at positive and 
negative angles Then, 7 \V/ l since includes 
j () The sums of th irst two even powers ol 


The upper limiting value, / of the kurtosis index 


for a complete diffusor as .V/ increases without bound 


The lower limit occurs when \/ that is when 
observations are taken at the specular angle and at 
| In this case ag 1LS—O3 Ld Thus we may 
state, mathematically 
lo a, S 
For the number ol the interva's ised in this 


that is five intervals on either side of the 
1] V/ 7 Thus for i completely 


that Is, On having 


analy SIs 
spec ular angle / 


diffusing surface constant 7 


a,=1.80— =~ 73 


Should we wish to define a seale of pe akedness for 
which a perfect diffuser is assigned the value zero 
then we may subtract a, from as and call the result 


ing difference the pe akedness index, 7’ 


Thus 


If the ordinates of one distribution are a constant 
factor of the corresponding ordinates of another 
distribution, then the indexes ¢, a,, and a, for one 
distribution will equal the corresponding indexes for 
the othe To distinguish one curve from the other. 
one point on each curve needs to be specified A 
convenient point is the value of the fractional 
reflectance at the specular ang 

The dispersion index, ¢. is at limit for 
completely spec ular reflection and 1s large for com- 
pletely diffuse reflection. The lower limit is zero, 
at the specular angle. The value 
diffuse determined as 





e 


its lowe 


because d 0 
for completely 
follows 


reflection is 





Vf SM —1\(/2M 
i; 


2\f+1 W(M+1 2M 
. lo 





the first \/ natural numbers may be computed fr 
the following equations: 


i as 2M+1 WiM-+1 
and 

wa 3477+3M—1 2M+1 W(M-+1 
3x (FA )*(FA (AS) 
Thus 

| VWiM 

i) ‘ 
l WM! 5 5M(M+1 


, i came 


For completely diffuse reflection, / constant 


Sut from eq 11 
255,.X ) 
n 


and from eq 12 


) 1 


Here also 7 2M+-1, 


( 


WiM+1 


» 


For the number of the intervals used in this analysis 
VW=5. Thus when / 
The characteristics for the indexes, 
to conclusions about their utility in describing su 
face appearance. Consider the ranges of these two 
indexes for the number of observations used in this 
Mathematically stated, from flat to 
ranges are 1.78<a,<@ and 


is constant, ¢=—3.16 
a, and a, lea 


analysis. 
peaked curves, the 
3.16>¢>0 

When a reflection curve is flat, a, has a small value 
in its range, while o has a large value in its rang 
Conversely, when a a completely 
specular reflection, a, is at its upper limit, and ¢ is 
at its lower limit. Thus a, may indicate glossiness 
directly, and o may indicate glossiness inversely. 


curve describes 


5. Application to Glossimetry 


The practice of glossimetry is now quite w 
established and generally vields satisfactory ais 
criminations, but must be recognized as an abridge 
goniophotometric technique. Just as in the pro! 


lems of colorimetry, the standardization of abridge: 
methods of measurement is necessary to enable o1 


to obtain on a common basis rapid compariso! 
The get 


during production of a specific material. 
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permitted to 


method of analyzing goniophotometric distribu- 

irves described In this paper may be the means 
ef arriving at instrumental geometry and 
recommended procedures for making objective evalu- 


some 


al of the surface appearance attribute ascribable 
~ mective observations on the specular gloss ol 
mens are usually made without use of angle- 
easuring devices Therefore, as the skewness 
ndex is a measure of dissymmetry of a distribution 
no pt is made to correlate this index with the 


ippearance attribute ascribable to gloss 
first problem confronting this analysis is, do 


ndexes a, and o correlate with surface appear- 
and distinctness-of-image gloss 
than the existing glossimetry 


the correlations prove satistactory 


ilar gloss, 
vell as ol bette 
chniques?§ If 
rroblem is how to int rpret this analvsis so 


sper 


he? x 
that iment geometry for abridged methods may 
ed wisely Three important geometric vari- 
a glossmeter are angle of incidence, (2 
aperture siz }) the ratio of the angular 
of the source to the receiver aperture 


5.1. Correlation with Distinctness of Reflected Images 


One of the attributes of surface appearance Is the 
stinctness with which images are reflected This 
mute s influenced by at least two properties of 


; mens. One of the properties is the smoothness 
of tl surface, the other is the diffuse refiection of 
the men Lack of smoothness tends to distort 

ted images, and diffuse reflection controls 


oht and shadow 
method has 
distinctness-of-image 
16] in 1936 described a photographic 
reflected 
was 


rast betwee n highlhi 


HH tolore no satistactory been de- 
to measure objective ly 
Hunter 
ehnique lor 


recording distinctness of 


However, no analysis 


ed in his paper. 


quantitative 


a. Experimental Conditions 


\s re is need for a quantitative index of dis- 
tinct ss ol reflec ted images the kurtosis index was 


tested { correlation with distinctness of image 
estin . For this experiment estimates were made 
by six observers, ranging in familiarity with the 


problems from thoroughly experienced to completely 
halve 
The 22 specimens used 
ials: Polished black 
ass, ae id-etched blac k glass polished white Vitro- 
white Vitrolite, and white cerami 
specimens placed on a black 
table. The table was placed near a 
ndow so that the specimens were illuminated by 
h-sky light. The made bv the ob- 
rs were on the distinctness of the unages of 
crossbars of the window Observers were 
turn, tilt, or 
ng the evaluations. The observers arranged the 
cimens in the order of estimated decreasing 
stinctness of the reflected images. The resulting 


following 
sand-blasted black 


were of the 
nter class, 
depolished 

These 


oth-cove red 


were 


estimates 


raise the specimens 





data are shown in table 3, together with the ranking 
obtained by the 


and 60 


kurtosis index, dispersion index, 


specular-gloss values of the specimens 


b. Rank Correlation 


Correlation indexes wert computed ou the basis 
ol Spearman's rank correlation coefficient p This 
correlation coeflicient has been found quite useful 
in correlating ordered o1 ranked data rhe equa- 


tion for computing p is 


where D is the rank number 
of items 


This 


difference, and N is the 


whether the 


When- 


equation indicates by its sign 


correlation is direct ol 


Inverse 
ever there is a tie in rank, the positions are split 
among the different items in the ti 
Tat 5 Re » 7" s 
MMB DBJ HKHGWH IN. 
\ 
\ 
ii 
A4 ‘ : ‘ ‘ ; 


The data indicate that confusion exists among the 
estimates where the specimens have nearly the same 
but one is black and the othe 
made 


and 


distinctness of image 
white. Therefore, 
for the entire group and for the group of blac 
the group of white specimens 

The correlation indexes by the SLX 
observers for the three correlation conditions are 
shown in table 4 These results indicate that 
relation of image distin tness with a, is better than 
and is ‘far better than with 60 specular 
When the specimens are separated into two 
groups according to reflectance, then correlation 
with a, is nearly perfect We may conclude that 
the diffuse reflectance of a specimen is an important 
factor in distinctness-of-image estimates 


computation Of p Was 


] 
obtained 
Core 


W ith 0 


gloss. 
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5.2. Correlation With 60° Specular Gloss 


Qne of the satisfactory glossimet 


the ht} specular orl ss method 


tained by analvsis of goniophotome 
relate with 60 specular gloss and also vield more 
information, then this analysis 
over the abridged method To 
exists between kurtosis and 60 


plot a, ona logarithm scal 


gloss readings for the 22 specimens listed in table ; 
fit the resulting 
points Figure »sa plot of the data show ing divi- 


and determine if a smooth curve 


sion mto two distinet groups shown 
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mimportan 
th i < 


we ad on of the 








depending on the appearances of the reflected imag 
Those specimens forming sharp images shaded 


Varving haze fall into one group The se specime 


lorming no sharp image s but having Varving spe 
lar gloss, fall in the other group. This examinat 
for correlation indicates not only that direct corre 


tion exists, but also that separation into disti 


groups in accord with visual appearance is acco 


plished by the coniophotometris eurve analy 
whereas the glossmeter readings alone do not suge 
the existence of these two groups 


Correlation of dispersion, o, with 60 specul 


gloss is shown in figure 4. This figure shows 
division of the specimens into the same two gro 
The correlation for o is inverse, that is, where 
oss is low, and conversely 





lar re, g 
We must conclude from the above data that 
analvsis not only correlates with existing glossims 
techniques but vields additional information abe 


surface appearance 


5.3. Selection of Instrumental Geometry 


This analysis has been shown lo have advanta 
over abridged methods in surface-appearance ev: 
ation Abridged methods of measurement us 
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vlossmeters, are required to obtain information econo- 
mically during the production of materials 
as well as goniophotometers, should have 


( 1loss- 


meters 

geometry to yield optimum discrimination between 
specimens. Assuming that the receiver is adequately 
sensitive, the resolution and the incidence angle may 


be chosen to vield this optimum discrimination 

['o determine the optimum geometric conditions 
for goniophotometry and glossimetry a of 
goniophotometric curves should be obtained with 
various materials for incidence angles ranging from 
perpendicular to grazing. Then the information 
from the goniophotometric curves needs to be cor- 
The data for 


series 


elated with visual estimates of gloss 
two specimens, A (acid-etched black glass) and S 
sand-blasted black glass), having 60° specular-gloss 
values of about 80, will serve to illustrate the pro- 


cedure for finding optimum geometry (see fig. 5). 
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f 
4 -2 re) 2 4 —4 ~2 1°) 2 4 
d, DEGREES 
Fic RE 5 ( om parison of 60 with 20 incidence ingles. 
60° incidence on left; 20° incidence on right; sand-blasted black glass, S; acid 


etched black glass, 


a. Choice of Incidence Angle 


Let us consider first how the peakedness index, P, 
may be used in the problem of choice of incidence 
angle 

Greatest discrimination between specimens is in- 
dicated where the difference between peakedness in- 
dexes is greatest. In the illustration, for 20° inci- 
dence, the difference between Ps (index for specimen 
S) and P, (index for specimen A) is Ps—P,=153 


6=151.4. For 60° incidence this difference is 
P;s—P,=379—2.8=376.2. Thus for greater dis- 


crimination between S and A, 60° incidence is super- 
ior to 20° incidence for goniophotometric measure- 
ments. 

This reasoning may be applied to glossmeter 
measurements. The receiver aperture of a gloss- 
meter includes flux reflected near the specular angle. 
Thus for optimum discrimination between speci- 
mens S and A, glossmeter measurements should be 








made at 60° rather than at 20 
that 


incidence, provided 
the flux accepted for measurement is suffi 


ciently close to the specular angle. To answer the 
question, what is meant by “sufficiently close’, we 
must consider the problem of choice of receiver 


aperture size. 
b. Choice of Receiver Aperture Size 


Let us now consider how the goniophotometric 
curves may be used in the problem of choice of re- 
ceiver aperture size. 

The goniophotometric curves provide sufficient 
information for the choice of receiver aperture size 
needed for discriminations. For goniophotometri: 
measurements the receiver aperture 
sufficiently small to trace adequately the geometric 
distribution of the fractional reflectances 
meters, on the other hand, integrate the flux enter- 
ing the receiver entrance window at the specular 
angle. The glossmeter must distinguish between 
different specimens by integrating between the limits 
of the receiver entrance window the areas under the 
respective goniophotometric curves. The limits of 
integration may be evaluated from either of the fol- 
lowing relationships: 


needs to be 


Gloss 


"a 


F,id=K, (16) 


"a 

I 6d 
e a e a 
or 


* red=K* F.8d, 17) 
a ‘.. 


- 


where a is the limit of integration, F; and F, are 
fractional reflectances of two specimens, d is the 
departure from the specular angle, and A is an ex- 
perimentally determined constant. Equation (16) 
states that the difference between the integrated 
areas is constant. Equation (17 that the 
ratio of the integrated areas is constant. The con- 
stant A should be evaluated for specimens of similar 
materials from visual estimates of either gloss dif- 
ferences or gloss ratios. No attempt should be made 
to evaluate A for radically dissimilar specimen types 


states 


6. Conclusions 


Since the surface appearance attribute ascribable 





to gloss may be evaluated most advantageously by 
monoplane goniophotometric curves, a technique for 
analyzing these curves needs to be adopted The 
variables of classical statistics, dispersion, skewness, 
and kurtosis, have been studied in this connection, 
| and it is found that the kurtosis index in particular 
| summarizes in a meaningful way the information 
given by the goniophotometric curve. Correlation 
between 60° specular gloss and kurtosis index is 
fairly good, and indicates that the two measure 
similar properties of surfaces such as those studied 
Furthermore, correlation between distinctness of 
image and kurtosis index is a decided improvement 
over correlation between distinctness of image and 
60° specular gloss. Kurtosis index is therefore pro- 
posed for adoption as a standard technique for re- 
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ducing goniophotometric data to meaningful num- 


bers 


This index has been applied successfully for 


choosing instrumental geor. ctry for abridged ronio- 
photometry, that is, glossimetry 


J.H 
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